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FOREWORD 


This  report  contains  the  description  of  a propeller  acoustics  test 
facility  utilized  to  conduct  a number  of  propeller  performance  and 
acoustic  tests  over  the  period  from  1970  to  1976.  The  facility  was 
developed  by  the  Turbine  Engine  Division  and  the  Technical  Facilities 
Division  of  the  Air  Force  Aero-Propulsion  Laboratory,  Wright-Patterson 
AFB,  Ohio, under  Project  3066,  Task  12,  Work  Unit  02.  The  effort  was 
conducted  by  Paul  A.  Shahady,  AFAPL/TBC,  and  Sigmund  W.  Kizirnis,  AFAPL/TFE 
during  the  period  of  1 January  1970  to  31  October  1976.  2/Lt  Robert  M. 
McGregor  provided  major  computer  support  to  the  project. 

Special  appreciation  is  extended  to  Raymond  Allen  for  providing 
outstanding  instrumentation  support  to  the  facility.  Appreciation  is  also 
given  to  George  Medisch,  foreman  of  the  facility  test  crew  during  most 
of  the  testing  period,  to  Wellington  Steel,  William  Howerton,  Harold  Crouch, 
Walt  Stebel,  and  Harold  Lee,  members  of  the  facility  test  crew  and  to 
Joseph  Simmons  for  his  support  to  the  acoustic  tests. 


TABLE  OF  CONTENTS 


SECTION 

I INTRODUCTION 

II  FACILITY  DESCRIPTION  AND  PERFORMANCE  MEASUREMENT 

SYSTEM 

III  ACOUSTIC  DATA  ACQUISITION  SYSTEM 

IV  ACOUSTIC  DATA  ANALYSIS  SYSTEM 

V ACOUSTIC  DATA  PRESENTATION  SYSTEM 

VI  FACILITY  ACOUSTIC  CHARACTERISTICS 

VII  CONCLUSIONS  AND  RECOMMENDATIONS 

VIII  REFERENCES 

IX  APPENDICES 

A.  DETAILED  SCHEMATICS  - ACOUSTIC  DATA  SYSTEM 

B.  DATA  PRESENTATION  COMPUTER  PROGRAMS 


LIST  OF  ILLUSTRATIONS 


FIGURE  page 

1.  Propeller  Acoustics  Test  Facility  at  Wrlght-Patterson  AFB  3 

2.  Performance  Data  Sheet  4 

3.  Typical  Microphone  Array  for  Propeller  Acoustics  Testing  8 

4.  Data  Flow  Schematic  for  Acoustic  Data  Acquisition  9 

5.  Data  Flow  Schematic  for  On-Line  Acoustic  Data  Analysis  11 

6.  Data  Flow  Schematic  for  Off-Line  Acoustic  Data  Analysis  12 

7.  Dynamic  Data  Analysis  System  at  the  Flight  Dynamics  Laboratory  13 

8.  Computerized  Performance  Input  Data  15 

8a.  Basic  Digitized  1/3-Octave  SPL  Data  16 

9.  Room  Resonances  in  Third  Octave  Bands  Around  Octave  Band 

Center  Frequencies  21 

10.  Uncertainty  Range  for  Radiation  from  a Simple  Sound  Source 

in  a Reverberant  Field  22 

11.  Relative  Sound  Pressure  Level  as  a Function  of  Distance  for  a 

Nondirectlonal  Point  Source  (Bldg  2QA)  24 

12.  Building  20A  Ambient  Noise  Levels  26 

13.  Rig  #2  Background  Noise  for  Drive  Motor  and  Propeller  Hub 

Without  Blades  27 

14.  Geometry  for  a Propeller  as  Viewed  from  a Measurement  Point, P, 

on  the  Axis  29 

15.  Sound  Level  Distribution  Vs  Distance  from  a Circular  Plane 

Source  31 

16.  Assumed  Directivity  Pattern  for  Sound  Radiation  from  Elemental 

Sources  Distributed  33 

17.  Theoretical  Noise  Patterns  for  Propeller  Noise  (Reference  9)  34 

18.  Geometry  for  an  Assumed  Elemental  Source  Directivity  Pattern 

as  Viewed  from  a Measurement  Point,  P,  on  the  Axis  35 

vl 


LIST  OF  ILLUSTRATIONS  (CWT 'D) 


FIGURE  PAGE 

19.  Sound  Level  Distributions  Vs  Distance  from  Circular  Plane 

Sources  With  Assumed  Directivity  Patterns  37 

20.  Relative  Sound  Pressure  Level  Vs  a Function  of  On-Axis  Distance 

from  a Propeller  of  Radius  R Meters  39 

21.  Geometry  for  a Propeller  as  Viewed  from  a Measurement  Point,  P, 

in  the  Plane  of  Rotation  40 

22.  Sound  Level  Distribution  from  a Circular  Arc  Source  of  Radius,  R, 

Vs  Distance  from  the  Center  in  the  Plane  of  the  Arc  43 

23.  Relative  Sound  Pressure  Level  for  a Propeller  of  Radius,  R Meters, 

Vs  Distance  from  the  Center  in  the  Plane  of  Rotation  44 


vli 


I.  INTRODUCTION 


The  purpose  of  this  report  is  to  describe  the  development  of  a 
research  facility  to  investigate  fan  and  propeller  noise  generation 
mechanisms.  Over  the  past  several  decades,  a considerable  amount  of 
propeller  noise  experimental  data  has  been  obtained.  Much  of  this 
data  was  acquired  using  simplified  test  rigs  employing  reciprocating 
or  gas  turbine  engines  as  the  means  for  driving  the  propeller.  Noise 
from  these  drive  systems  seriously  contaminated  much  of  the  measured 
propeller  noise  under  certain  operating  conditions.  The  more  recently 
constructed  turbine  engine  fan  test  facilities  have  alleviated  the 
drive  noise  contamination  problem  by  properly  suppressing  drive  system 
noise.  However,  most  of  the  facilities  capable  of  large-scale  acoustic 
testing  are  located  outdoors  and  are  subject  to  weather  conditions. 
Therefore,  a program  was  initiated  to  develop  an  all-weather  test 
capability  to  measure  near-  and  far-field  propeller  and  fan  noise 
free  of  contamination  from  excessive  background  noise. 

A 3500- horsepower  electric  whirl  rig,  located  in  Building  20A 
at  Wright-Patterson  AFB,  was  modified  to  allow  both  near-  and  far- 
field  acoustic  testing.  The  rig  is  located  in  a large  building  with 
acoustically  treated  walls.  Reverberation  time  and  sound  pressure  level 
measurements  were  conducted  to  determine  the  acoustic  characteristics 
of  the  building.  Acoustic  data  was  acquired  using  a 12-channel  portable 
data  acquisition  system  consisting  of  microphone  stands,  microphones, 
land  lines,  signal  conditioning  equipment,  a multi-channel  analog  FM 
tape  recorder,  and  associated  calibration  equipment.  Acoustic  data 
is  frequency  analyzed  into  one- third  octave  or  narrow  bands  and  auto- 
matically plotted  using  a computerized  dynamic  analysis  system  located 
in  the  Air  Force  Flight  Dynamics  Laboratory  at  Wright-Patterson  Air 
Force  Base.  A series  of  data  presentation  programs  have  been  developed 
to  cross  plot  a number  of  key  performance  and  acoustic  parameters. 


II.  FACILITY  DESCRIPTION  AND  PERFORMANCE  MEASUREMENT  SYSTEM 


The  3500-horsepower  electric  whirl  rig  consists  of  a large 
concrete  pier  which  houses  the  electric  drive  motor,  thrust  and  RPM 
measuring  equipment,  and  various  accessory  drives.  The  pier  rises 
about  25  feet  off  the  floor  of  a large  open  building.  The  facility 
side  walls  and  floor  were  covered  with  six-inch  thick  Coustic  TM 
polyurethane  foam  to  minimize  acoustic  reflections.  The  control 
room  is  located  under  the  pier  and  the  propeller  may  be  observed 
via  a periscope  and  strobe  light  arrangement.  Closed  circuit  TV 
cameras  are  also  used  to  monitor  the  facility.  Various  parts  of  the 
facility  are  shown  in  Figure  1.  The  following  technical  summary  of 
the  performance  measurement  system  is  taken  from  Reference  1.  The 
input  power  to  the  propeller  is  calculated  from  measuring  the  armature 
voltage  and  amperage  at  the  electric  drive  motor.  Predetermined 
correction  factors  are  then  applied  to  allow  for  the  copper  and 
field  winding  losses.  The  resultant  watts  are  then  converted  to 
horsepower  and  an  atmospheric  correction  factor  is  used  to  adjust 
the  data  to  standard  day  conditions.  These  calculations  are  made 
with  an  ei  >'onic  desk  calculator  during  the  course  of  the  test  so 
that  fip  ves  can  be  drawn  to  check  for  obvious  data  discrepancies. 
The  sy  '■alibrated  against  a well  documented  test  propeller 

whose  e been  accurately  set  and  locked  into  place.  No- 

load  are  determined  by  motoring  the  rig  at  various  RPM 

settings  without  a propeller  attached. 

The  thrust  is  measured  by  converting  the  movement  of  the 
propeller  shaft  to  hydraulic  pressure  via  a hydraulic  diaphragm. 

The  pressure  signal  is  then  directly  converted  to  pounds  thrust  with 
a precalibrated  Emery-Tate  load  indicator  and  then  corrected  to 
standard  day  conditions.  The  thrust  system  is  calibrated  statically 
by  applying  a known  load  (lead  weights)  to  the  propeller  shaft. 

Accurate  shaft  RPM  was  obtained  from  a magnetic  pickup  which  receives 
impulses  from  the  drive  motor  shaft.  These  impulses  are  then  presented 
on  a digital  display  in  the  control  room  as  propeller  RPM. 

All  of  the  previously  mentioned  performance  data,  horsepower, 
thrust,  and  RPM,  is  taken  directly  from  the  instrumentation  in  the 
control  room.  This  data  is  then  corrected  for  the  various  loss  and 
atmospheric  effects  by  employing  pre-determined  correction  factors. 

The  data  is  tabulated  on  the  form  shown  in  Figure  2 and  turned  over 
to  the  project  engineer  for  reduction.  This  is  accomplished  by  using 
a computer  program  entitled  "Computer  Program  for  Reducing  Static 
Propeller  Test  Data"  (Reference  2).  The  following  explanation  of 
the  program  generally  follows  that  of  Chopin  (Reference  3). 

The  program  accepts  whirl  rig  data  in  the  format  in  which  it  is 
taken  from  the  rig  by  the  test  crew.  This  is  then  reduced  by  the 
computer  into  power  coefficients  (C  ),  thrust  coefficients  (C^), 
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PROPELLER  WHIRL  TEST  DATA  SHEET 


Figure  2:  Performance  Data  Sheet 


of  merit  (F.M.),  thrust/horsepower  (Th/HP),  and  propeller  tip 
The  output  format  of  the  program  is  shown  in  Table  1.  The 
symbols  presented  in  Table  1 are  defined  as  follows: 


C /C  , figure 
Mach^number, 


PROPELLER  CHARACTERISTICS 
BETA  - Test  blade  angle 

AF  - Blade  activity  factor 

DIA  - Propeller  diameter  in  feet 

NBL  - Number  of  blades 

TEMPC  - Ambient  temperature  in  degrees  Centigrade 

TEMPR  - Ambient  temperature  in  degrees  Rankine 

SIGMA  - Density  ratio 


RAW  DATA  POINTS 

RPM  - Propeller  rpm 

HP  - Corrected  horsepower 

TH  - Corrected  thrust 

TMACH  - Propeller  tip  Mach  number 
RCT  - Raw  thrust  coefficient 

RCP  - Raw  power  coefficient 

RCT/CP  - Ratio  of  raw  thrust  to  raw  power  coefficient 

RFM  - Raw  figure  of  merit 

RTH/HP  - Ratio  of  raw  corrected  thrust  to  corrected  horsepower 


FITTED  CURVE 
MACH 
HP 
TH 

TIPS 

RPM 

CT 

CP 

CT/CP  - 
RM 

TH/HP  - 


DATA  FOR  CONSTANT  MACH  NUMBER  INCREMENTS 
Selected  Mach  number  increment 
Horsepower  at  Mach  increment 
Thrust  at  Mach  increment 

Propeller  tip  speed  in  ft/sec  corresponding  to  Mach  increment 

Propeller  rpm  at  Mach  increment 

Thrust  coefficient  at  Mach  increment 

Power  coefficient  at  Mach  increment 

Ratio  of  thrust  coefficient  to  power  coefficient  at  Mach 
i ncrement 

Figure  of  merit  at  Mach  increment 

Ratio  of  thrust  to  horsepower  at  Mach  increment 


TABLE  1:  PROPELLER  PERFORMANCE  DATA  OUTPUT  FORMAT 
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III.  ACOUSTIC  DATA  ACQUISITION  SYSTEM 


The  acoustic  data  acquisition  system  consists  of  Bruel  and  Kjaer  half- 
inch microphones  and  cathode  followers  and  in-house  designed  line  drive 
amplifiers  mounted  on  portable  microphone  stands.  The  microphones  are 
located  23  feet  above  the  floor  in  a horizontal  plane  through  the  propeller 
centerline.  A typical  microphone  array  used  for  several  test  programs  is 
shown  in  Figure  3. 

The  data  signals  are  carried  by  250-foot  land  lines  to  an  instrumentation 
control  room  located  below  the  rig.  The  signals  were  conditioned  through 
variable  gain  amplifiers  and  recorded  on  a 14-channel  AMPEX  FR  1300  tape 
unit  for  later  analysis.  A data  flow  schematic  of  the  acoustic  data 
acquisition  system  is  shown  in  Figure  4.  The  two  symbols  for  the  tape 
recorder  represent  the  data  flow  in  the  record  mode  and  in  the  reproduce 
mode  respectively.  Gain  and  attenuator  adjustments  are  accessible  at 
various  points  in  the  system.  Phase  and  gain  controls  are  located  at  the  line 
drive  amplifier  for  system  calibration  adjustments.  A zero  to  twenty  dB 
attenuation  is  available  at  the  junction  box.  Final  adjustments  are  made  at 
the  variable  gain  amplifiers  before  the  signals  are  input  to  the  tape  recorder. 
The  data  signals  can  be  monitored  at  various  points  in  the  system.  A signal 
output  BNC  is  located  at  the  line  drive  amplifier  and  at  the  junction  box. 
Individual  channel  scopes  are  used  to  monitor  either  the  tape  recorder  input 
signals  or  the  reproduced  output  signals.  A signal  selector  circuit  is  used 
to  provide  any  of  the  tape  recorder  input  signals  for  detailed  analysis.  These 
signals  are  recorded  on  the  loop  recorder  for  on-line  frequency  analysis.  These 
various  monitor  outputs  allow  the  signal  to  be  evaluated  on  the  output  side  of 
each  major  component  from  the  line  drive  amplifier  to  the  tape  recorder.  The 
total  system  is  calibrated  at  the  beginning  and  end  of  each  tape  and  at  the 
beginning  and  end  of  each  test  day.  An  NBS  traceable  piston  phone  provides 
a 124  dB  pressure  signal  at  250  Hertz  to  the  microphone  diaphragm.  This 
signal  is  recorded  on  tape  for  each  microphone  channel.  Gains  for  the 
various  data  acquisition  system  components  are  recorded  manually  on  data 
log  sheets. 

The  data  acquisition  system  is  primarily  comprised  of  off-the-shelf 
general  purpose  instrumentation.  However,  four  of  the  components,  the 
line  drive  amplifier,  junction  box,  signal  selector,  and  audio  monitor 
system  were  designed  in-house.  Detailed  schematics  of  these  components  are 
presented  in  Appendix  A. 
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FIGURE  3: 

TYPICAL  MICROPHONE  ARRAY  FOR 
PROPELLER  ACOUSTICS  TESTING 


FIGURE  4:  DATA  FLOW  SCHEMATIC  FOR  ACOUSTIC  DATA  ACQUISITION 


IV.  ACOUSTIC  DATA  ANALYSIS  SYSTEM 


There  are  several  systems  available  for  analysis  of  acoustic  data. 

One  system  available  for  single  channel  on-site  analysis  consists  of  an 
AMPEX  FR  1300  multichannel  tape  recorder,  scopes,  variable  gain  amplifiers, 
an  AMPEX  loop  recorder,  and  B&K  frequency  analysis  instrumentation.  A data 
flow  schematic  of  this  acoustic  analysis  system  is  shown  in  Figure  5. 

Signals  are  reproduced  from  the  tape  recorder  and  are  monitored  on  individual 
channel  scopes.  They  are  then  input  to  variable  gain  amplifiers.  Output 
signals  from  the  amplifiers  are  routed  either  to  the  loop  recorder  and  back 
to  a signal  selector  or  to  the  signal  selector  directly.  The  signal  is  then  a- 
vailable  at  the  output  of  the  selector  circuit  for  either  narrow  band  or  one- 
third  octave  analysis.  The  signals  can  be  monitored  aurally,  visually,  and 
graphically  using  the  audio  monitor,  an  oscilloscope,  and  a level  recorder. 

A second  system  is  available  off-site  at  the  Air  Force  Flight  Dynamics 
Laboratory.  The  system  handles  third-octave  data  in  real  time  using  a 
General  Radio  Model  1933  real  time  1/3  octave  analyzer.  Results  are  available  in 
tabular  output  or  in  plot  tape  form.  Narrow  band  analysis  is  accomplished  using 
a digital  FFT  system.  The  data  flow  schematic  for  both  systems  are  presented 
in  Figure  6.  The  analog  signal  from  the  tape  deck  is  first  edited  and  then 
digitized  in  frequency  decades.  Narrow  band  analysis  is  then  conducted  using 
a digital  FFT  system.  The  results  are  tabulated,  run  through  a sort/merge 
process  and  a post  processor  to  create  a plot  tape.  The  results  are  then 
plotted  on  a calcomp  plotter.  Pictures  of  the  Air  Force  Flight  Dynamics 
Laboratory's  data  analysis  system  are  shown  in  Figure  7. 


FIGURE  5:  DATA  FLOW  SCHEMATIC  FOR  ON-LINE  ACOUSTIC  DATA  ANALYSIS 


A 
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Figure  7:  Dynamic  Data  Analysis  System  at  the  Flight  Dynamics  Laboratory 


Several  computer  programs  were  written  to  reduce  and  plot  propeller 
performance  and  acoustic  test  data.  The  data  comes  in  two  basic  forms: 
Propeller  Whirl  Test  Data  Sheets  (AFAPL  Form  14A),  and  digitized  sound 
spectra  computer  cards. 

Data  from  the  test  rig  sheets  must  first  be  coded  onto  computer  cards 
as  shown  in  Figure  8.  In  this  format  it  is  compatible  with  programs  written 
for  data  reduction  programs  in  ASD-TR-68-19,  "Computer  Program  for  Reducing 
Static  Propeller  Test  Data."  Not  all  data  on  these  cards  will  be  required 
by  each  individual  program,  but  with  commonality,  the  number  of  cards  and 
the  time  spent  punching  them  is  reduced. 

The  digitized  spectra  data  cards  were  obtained  from  a General  Radio 
1933  Real  Time  third-octave  band  analyzer,  and  come  in  two  formats.  The 
first  format  covers  the  frequencies  from  3.15  Hz  through  20,000  Hz,  while 
the  second  covers  20  Hz  through  20,000  Hz.  Each  block  of  data  consists 
of  four  cards,  the  first  containing  the  run  and  microphone  channel 
identification  numbers.  The  last  number  on  the  last  card  is  the  overall 
sound  pressure  level  (OASPL)  in  dB  for  that  block.  All  other  numbers 
represent  the  SPL's  for  their  respective  third-octave  bands.  Figure  8. a. 
shows  data  sets  for  three  microphones  (4,  5,  6)  during  run  number  six  of 
a particular  configuration. 

The  basic  coded  rig  data  and  spectra  data  are  used  by  five  programs: 

HORSE,  THORSE,  NOISYl,  N0ISY2,  and  DBA.  Each  program  is  written  in  FORTRAN 
IV-Extended  for  the  CDC  6600  computer,  and  with  the  exception  of  DBA,  requires 
a CALCOMP  plotter  and  associated  software.  DBA  requires  a peripheral  card 
punch.  References  2 and  4 were  used  in  the  development  of  these  programs. 


HORSE 

Program  HORSE  plots  thrust  vs.  horsepower  for  any  configuration.  Several 
parameters  which  control  the  output  format  have  been  preset  within  the 
main  program  body.  Changes  should  be  made  using  a NAMELIST  card.  Input 
may  include  600,  750,  900,  1050,  and  1200  RPM  data  or  600,  900,  and  1200 
RPM  data.  If  five  points  are  entered,  all  five  will  be  used.  Likewise, 
if  three  are  entered,  three  will  be  plotted.  Setting  the  Namelist  variable 
lORPM  = 2 gives  three-point  output  even  if  five  points  are  entered.  Three 
blade  angles  are  required  for  each  RPM.  Spline-fitted  curves  through 
all  points  are  standard.  Thrust  and  horsepower  limits  are  preset  at  1600 
and  240  respectively,  and  are  held  constant  for  easy  comparison  between 
plots.  If  desired,  the  user  may  specify  his  own  limits.  Different  limits 
need  only  be  set  once  for  any  number  of  plots  generated  during  a single 
program  run.  If  changes  are  required,  however,  99999  in  columns  1 - 5 
alerts  the  program  to  a new  $INPUT  card.  A program  listing  with  sample 
input  cards  and  output  plots  is  found  in  Appendix  B.1. 
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FIGURE  8a:  BASIC  DIGITIZED  1/3-OCTAVE  SPL  DATA 


THORSE 


Program  THORSE  is  a derivative  of  HORSE,  and  plots  the  ratio  of  thrust 
to  horsepower  versus  horsepower.  Input  cards  are  exactly  the  same  as  those 
for  HORSE.  Any  section  of  a plot  which  falls  below  the  X-axis  is  automatically 
eliminated.  This  is  not  true  for  values  exceeding  the  Y-axis  maximum  (preset 
at  11.)  where  the  user  risks  running  the  plotter  into  the  travel  stops,  and 
throwing  all  remaining  points  out  of  place.  A program  listing  with  sample 
input  cards  and  output  plots  is  found  in  Appendix  B.2. 


NOISYl 

Program  NOISYl  plots  SPL  (dB)  versus  third-octave  band  frequency.  The 
user  may  select  any  range  between  20  Hz  and  10,000  Hz  for  plotting.  The 
overall  SPL  for  the  selected  range  is  plotted  separately  on  the  right  of  the 
graph.  The  user  may  elect  to  have  the  A-weighted  spectra  plotted  along  with 
the  dB  spectra.  Any  number  of  spectra  may  be  put  on  the  same  axes  as  long 
as  they  are  from  the  same  microphone  channel.  As  in  the  previous  programs, 
plotting  parameters  are  preset  but  may  be  changed  using  the  standard  NAMELIST 
format.  A blank  card  at  the  end  of  a block  of  data  signals  an  end  to  plotting 
on  that  set  of  axes.  The  user  then  has  three  options:  (1)  end  that  session 

of  plotting  with  an  ENDDATA  card,  (2)  start  another  plot  with  a new  set  of 
data  cards  or,  (3)  change  the  output  Namelist  parameters  using  a NEWPARAM  card. 
Following  the  NEWPARAM  card  should  be  a new  set  of  data.  The  program  listing, 
sample  input  cards  and  output  plots  are  found  in  Appendix  B.3. 


N0ISY2 

Program  N0ISY2  is  a derivative  of  NOISYl.  N0ISY2  plots  the  spectra 
of  several  microphones  on  the  same  axes  as  long  as  the  RPM  remains  constant. 
This  differs  from  NOISYl  where  the  RPM  varied  for  a single  microphone. 

A program  listing  with  sample  input  cards  and  output  plots  is  found  in 
Appendix  B.4. 


I DBA 


Program  DBA  reduces  the  digitized  spectra  data  to  dB  and  dBA  sound 
• pressure  levels.  The  program  sorts  through  each  run  of  the  basic  data, 

f stopping  at  user-specified  microphone  numbers.  Over  any  given  spectra 

I range,  the  program  will  then  calculate  the  overall  SPL  in  both  dB  and 

dBA.  The  data  are  then  punched  on  cards  with  the  corresponding  identifica 
tion  numbers  and  calculation  range  indices.  If  cards  are  not  desired, 

' the  user  may  set  IPUNCH  = 0 and  receive  only  the  printput.  Punched  cards 

from  DBA  are  used  by  program  NEWCITY.  A program  listing  with  sample  input 
cards  and  an  output  listing  are  found  in  Appendix  B.5. 


NEWCITY 

The  final  data  presentation  program  is  NEWCITY.  Using  both  the 
basic  r1g  data  and  the  SPL  data  cards  from  DBA,  NEWCITY  will  plot  SPL  vs. 
thrust  and  SPL  vs.  horsepower.  Several  options  for  presentation  are  open 
to  the  user.  The  user  may  elect  to  receive  either  the  thrust  or  horsepower 
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plots  or  both.  The  legend  block  may  be  eliminated.  Three  or  five  RPM 
values  may  be  used,  depending  on  the  input  data.  Output  curves  may  be  in 
either  dB  or  dBA,  with  spline  or  up  to  6th-order  polynomial -fitted  curves. 

If  the  number  of  points  is  too  few  for  the  requested  polynomial  fit,  NEWCITY 
will  default  to  the  highest  possible  polynomial.  A program  listing  with 
sample  input  cards  and  an  output  listing  are  found  in  Appendix  B.6. 


VI.  FACILITY  ACOUSTIC  CHARACTERISTICS 


Ideally,  we  would  like  to  be  able  to  measure  noise  characteristics  in  a 
free- field  environment  uncontaminated  by  acoustic  reflections,  room  mode 
interactions,  reverberation,  background  noise,  etc.  Unfortunately,  we  are 
faced  with  all  these  problems  when  we  attempt  to  evaluate  acoustic  sources 
in  an  enclosure.  We  must,  therefore,  either  minimize  these  problems  or 
carefully  define  their  effect  on  the  acoustic  characteristics  of  the  source. 


Acoustic  Reflections 

The  problem  of  acoustic  reflections  contaminating  data  measured  on 
propeller  whirl  rig  #2  has  been  minimized  by  the  use  of  high  absorption 
acoustic  treatment  on  the  near  walls  and  the  floor  of  the  test  rig.  Six- 
inch  thick  pads  of  Custifoam  were  used  to  eliminate  all  reflections  above 
five  hundred  Hertz  and  minimize  reflections  in  the  100  to  500  Hertz  range. 


I 

i 

i 


Room  Mode  Interactions 

The  acoustic  power  output  of  a source  depends  on  the  impedance  presented 
to  it  by  the  surrounding  medium.  Reverberant  room  modes  may  increase  or 
decrease  this  impedance  depending  on  the  location  of  the  source  and  on  its 
frequency  relative  to  the  room  mode  frequencies.  If  the  room  is  very  large 
in  comparison  with  the  acoustic  wavelength,  there  will  be  so  many  modes  at 
any  given  frequency  that  the  effects  of  the  individual  modes  are  likely  to 
cancel  out.  The  total  number  of  room  resonances  occurring  in  a rectangular 
room  in  the  frequency  range  from  0 to  f is  given  by  the  expression. 


Q = 

Where: 


In  a narrow  frequency 
are 


3 

V 5 room  volume  (m  ) 


(1) 


S i room  surface  (m^) 

L = room  edges,  ly  ly 
c i velocity  of  sound  (m/sec) 

band  ( af  = fp  - f ^ ) around  f the  number  of  resonances 


aQ  - Q»  - Qf 
^2  ^1 

For  Building  20A,  V = 69,859  m^,  S = 13,012  m^  and  L 


(2) 

170  m 

i 
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Figure  9 gives  the  number  of  room  resonances  for  the  propeller  test  facility 
for  one-third  octave  bands  at  octave  band  center  frequencies  from  63  Hertz 
to  8 KHz.  Figure  10  shows  the  uncertainty  of  the  power  output  at  a single 
arbitrary  source  position  for  a simple  source  for  both  pure  tones  and  bands 
of  noise.  This  shows  that  the  expected  error  can  be  neglected  if  the  room  is 
large  enough  assuming  that  the  source  is  more  than  a wavelength  from  a wall 
(Reference  5).  For  the  propeller  rigs,  the  near  walls  and  the  floor  ad- 
versely affect  the  data  below  60  Hertz, 


Reverberation 


The  reverberation  characteristics  of  the  facility  were  evaluated  using 
the  Norris-Eyring  Formula  (Reference  6)  where  the  reverberation  time  in 
seconds  is  given  by  the  expression 


T = 


S[-2.3  log  (1-ff)] 


3 

Where:  V = volume  of  room,  m 

2 

S = area  of  bounding  surface,  m 
” = average  absorption  coefficient  of  the  room 
The  average  absorption  coefficient  is  therefore  given  by 


^ = 1 - 10-°'''/^^ 

The  room  constant,  R^,  is  given  by  the  expression 


S“ 

l-« 


(4) 


(5) 


Reverberation  time  measurements  were  made  in  the  facility  using  a rifle 
as  the  sound  source,  A number  of  shots  were  recorded  on  magnetic  tape  at 
various  locations  in  the  facility.  The  data  was  then  analyzed  in  one-third 
octave  frequency  bands  at  octave  band  center  frequencies.  Plots  of  amplitude 
versus  time  for  these  bands  were  used  to  determine  the  reverberation  time. 

The  reverberation  time,  average  absorption  coefficient,  and  room  constant 
are  given  as  a function  of  third-octave  bands  at  octave  band  center  frequencies 
in  Table  2.  Figure  11  gives  the  relative  sound  pressure  level  as  a function 
of  distance  from  the  acoustic  center  of  a non-directional  source  for  Building 
20A.  Relative  sound  pressure  level  is  given  by  the  expression 


SPL-PWl  . 10  )o,,„  * I-) 


Airr 


(6) 
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NUMBER  OF  ROOM  RESONANCES, 
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FIGURE  10:  UNCERTAINTY  RANGE  FOR  RADIATION  FROM  A SIMPLE  SOUND  SOURCE  IN  A REVERBERANT  FIELD 


FIGURE  11-  RELATIVE  SOUND  PRESSURE  LEVEL  AS 
A FUNCTION  OF  DISTANCE  FOR  A NON- 
DIRECTIONAL  POINT  SOURCE  (Blda  20A) 
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DISTANCE  FROM  ACOUSTIC  CENTER  OF  NON-DIRECTIONAL  SOURCE,  D (meters) 


2 

Where:  is  the  room  constant  in  m 

r is  the  distance  in  m 

Q is  a directivity  factor  (equal  to  1 for  propeller 
test  rigs) 

Since  the  room  constants  for  Building  20A  are  the  same  for  all  frequency 
bands  between  63  Hertz  and  2 KHz,  this  range  is  represented  by  the  top 
curve  in  Figure  11.  The  straight  line  shows  the  variation  of  relative 
sound  pressure  level  with  distance  for  a free-field  environment  (no 
reflections,  reverberations,  etc.). 


Background  Noise 

Background  noise  during  testing  in  the  Building  20A  facility  is  not 
a significant  problem.  Figure  12  shows  a plot  of  background  noise  without 
the  facility  drive  motors  operating.  Figure  13  shows  the  background  noise 
as  a function  of  RPM  for  a rotating  propeller  hub  without  blades.  Neither 
case  results  in  significant  background  noise  when  compared  to  the  noise 
generated  by  a typical  propeller.  The  measurements  were  made  at  a typical 
microphone  position  approximately  3 meters  from  the  source.  During  each 
propeller  test  program,  ambient  and  drive  motor  background  levels  are 
recorded  as  a function  of  RPM  for  each  microphone  position.  In  this  way 
the  source  noise  levels  can  be  compared  directly  with  the  background  noise 
to  ensure  that  the  signal -to-noise  ratio  is  adequate.  Generally,  we 
require  that  the  signal  be  at  least  10  dB  above  the  background  noise  in 
each  third-octave  band  of  interest. 

We  have  now  evaluated  several  problem  areas  relative  to  acoustic  testing 
in  Building  20A--acoustic  reflections,  room  mode  interactions,  reverberation, 
and  background  noise.  One  other  area  should  be  investigated  to  properly 
understand  the  acoustic  data  obtained  during  propeller  testing.  Since 
acoustic  propagation  characteristics  with  distance  change  depending  on  the 
source,  we  should  know  what  to  expect  in  terms  of  noise  level  variation  with 
distance  from  the  propeller.  It  is  beyond  the  scope  of  this  report  to 
discuss  all  of  the  noise  generation  and  propagation  theories  applicable  to 
propellers.  An  excellent  review  of  propeller  noise  theory  is  contained 
in  Reference  7.  It  Is  sufficient  to  point  out  that  the  broadband  character- 
istics of  a propeller  can  best  be  evaluated  from  noise  measurements  made 
along  the  axis  of  the  propeller.  Along  the  axis  there  Is  no  Doppler  effect. 
Also,  the  acoustic  spectrum  measured  on  the  axis  Is  generated  entirely  by 
fluctuating  loads  on  the  blade  and  Is  simply  related  to  their  integrated 
spectrum.  Measurements  made  In  the  plane  of  rotation  Involve  predominantly 
rotational  noise  due  to  harmonics  of  the  blade  passage  frequency. 

Reference  8 provides  the  basis  for  a preliminary  analysis  of  noise 
propagation  with  distance  from  a propeller.  In  our  application,  we  use 
a slightly  modified  version  of  Maekawa's  result  for  circular  plane  noise 
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1200  RPM 


FIGURE  13:  RIG  #2  BACKGROUND  NOISE  FOR  DRIVE  MOTOR  AND  PROPELLER  HUB  WITHOUT  BLADES 


sources  in  order  to  predict  propagation  characteristics  along  the  axis 
of  a propeller.  For  measurements  made  in  thi  plane  of  rotation  we  have 
developed  a new  analysis  based  on  propagation  from  a curved  array  of  point 
sources. 


Noise  Propagation  Along  the  Propeller  Axis 

When  viewed  from  the  axis,  the  propeller  can  be  represented  as  a circle 
of  radius  R made  up  of  point  sources  distributed  continuously  on  the  surface 
and  radiating  noise  energy  spherically  in  random  phase  so  that  the  nature  of 
wave  motion  can  be  neglected.  Sound  intensity  at  a particular  measurement 
point  will  be  obtained  by  integrating  the  sound  energy  radiating  from  each 
point.  Consider  the  sound  intensity  at  a measurement  point,  P,  on  the 
propeller  axis  as  shown  in  Figure  14.  The  energy  density  at  P,  a distance, 
d,  from  the  propeller  center  is  given  by 


E = 


W dS 


Js 


i(d  +x^)c 

Where:  W is  the  sound  power  per  unit  area 

Equation  (7)  can  be  rewritten  as 


(7) 


tanede 


Evaluating  the  integral  from  0 to  « we  obtain 
E 


(8) 


^ [-In  cos  e IJ]  * ^ [-In 


(R2.d2)V2| 


^ {-In  [1  + (d/R)-2]-l/2| 
^ {in  [1  + (d/R)"^]} 


We  can  express  this  in  terms  of  decibels  as  follows: 
E 


Where:  p is  the  sound  pressure  (N/m^) 

p is  the  density  of  air  (kg/m^) 
c Is  the  speed  of  sound  (m/sec) 


(9) 

(10) 
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(12) 


= W |ln  [1  + (d/R)'^]} 

2 -12 

By  referencing  the  pressure  to  .00002  N/m  and  the  power  to  10  watts 
we  obtain 


^ref 

Where: 


ref 


-12 

(f)  l— 7^  {in  [1  + (d/R)-2]l  (13) 

^ (.00002)'^  ' ' 


pc  = 408  kg/m^/sec  for  air 


Taking  10  log^g  of  each  side  we  obtain 

(14) 

2 

10  log^o  (^2 ) = 10  log^g  (g ) + 10  log^g  (.255)  + 10  {in  [1  + (d/R)"^] 

P pg^  10 

Equation  (14)  can  be  expressed  as 

SPL  = PWL^  -6+10  log^Q  {in  [1  + (d/R)‘^]|  (15) 

Where:  PWL^  is  the  power  level  per  unit  area 

Figure  15  is  a plot  of  the  third  expression  in  equation  (15)  as  a function  of 
d/R.  The  plot  shows  that  in  a free-field  environment,  the  noise  level  should 
drop  off  at  6dB  per  doubling  of  distance  for  all  points  along  the  axis  greater 
than  approximately  one  diameter  (2R)  away. 

Equation  (15)  was  developed  by  assuming  that  source  elements  on  the 
surface  of  the  propeller  radiate  omnidirectionally.  A similar  equation 
can  be  developed  by  considering  that  each  surface  element  on  the  plane  noise 
source  has  some  directivity  of  its  own.  A directional  radiation  power  can 
be  given  by  the  expression: 


cos"  6 (16) 

u 0 

Where:  J is  the  radiation  power  to  the  normal  of  a surface 

° element  given  in  the  units  watt/m^/sterad. 

e is  the  angle  between  Jg  and  the  normal  to  the  surface. 
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NORMALIZED  ON-AXIS  DISTANCE,  d/R 


Figure  16  shows  the  form  of  the  assumed  directivity  patterns  for  various 
values  of  n.  In  the  figure,  when  n = 0 the  element  is  omnidirectional  and 
when  n = « the  sound  is  radiated  only  to  the  normal  and  makes  a plane  wave. 
Figure  17  shows  the  characteristics  of  the  directivity  of  propeller  noise 
for  various  generation  mechanisms.  The  patterns  assumed  for  this  study  are 
representative  of  those  expected  for  on-axis  measurements. 


The  radiation  power  of  the  surface  per  unit  area,  W,  is  given  by 

j/2  ^Tf/2 


J-  2tt  sin6d0  = 2ttJ.  / 

0 ® “ ■'o 


cos  esinede  (17) 


Using  the  identity  /cos'^esinede  = -cos*^^^  /(n+1),  we  obtain 


u _ 2ttJ 
W = 0 


Using  equations  (7),  (8)  and  (18),  the  sound  energy  density  at  a 
receiving  point,  P,  as  defined  in  Figure  18  is  given  by 


= ~f 

ZcJn 


‘J„  2if  X sine  de 

b 


0 X cose 


= i f 


Jq  cos  0tan0d0 


When  n = 0,  equation  (20)  reduces  to 


= - /*  J 


tan0d0 


Noting  that 


J . 

0 2ir 


■ 5?  L 


tanedo 


the  same  as  equation  (8) 
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FIGURE  18; 


GEOMETRY  FOR  AN  ASSUMED  ELEMENTAL  SOURCE 
DIRECTIVITY  PATTERN  AS  VIEWED  FROM  A MEASUREMENT 
POINT,  P,  ON  THE  AXIS. 


When  n > 0,  equation  (20)  becomes 


^ rcos"o 


tano  do 


(24) 


Using  the  identity  fcos^e  tano  d©  = -cos'^e/n,  we  obtain 


E = ^ [1  - (cos.)"] 


(25) 


Replacing  cos“  with  the  proper  rectilinear  coordinates  yields 


^ = h ^ - (r  V’J  )"3 


V(R/d)‘ 


(26) 


Repeating  the  approach  used  to  derive  equations  (11)  through  (15),  we  obtain 


SPL  = PWL.  -3  + 10  log,.  - (~i (27) 

^ 10  / n V(R/d)  +1 


Figure  19  is  a plot  of  the  third  expression  in  equation  (27)  as  a 
function  of  d/R,  for  various  directivity  indices,  n.  The  directivity 
index,  n,  increases  as  the  frequency  of  the  source  element  increases.  There- 
fore, at  high  frequency,  sound  drop-off  with  distance  on  the  center  axis  is 
gradual.  This  is  not  true  for  measurement  points  out  of  the  center  axis. 

In  order  to  compare  equations  (15)  and  (27)  with  equation  (6), 
the  equations  must  be  converted  to  total  power  level  equations  where  the 
total  power  level,  PWL,  is  given  by 


PWL  = PWL^  + 10  log^Q  nR^ 


(28) 


Therefore,  equations  (15)  and  (27)  can  be  expressed  as  follows: 


For  n * 0 


SPL-PWL 


-6  - 10  log.|Q  ttR^  + 10  log.|Q 


ln[l+(d/R) 


(29) 
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NORMALIZED  ON-AXIS  DISTANCE,  d/R 


For  n > 0 


SPL-PWL  = -3-10  log^Q  ttR^  + 10  log^^  jXn^i-(  j (30) 


Figure  20  is  a plot  of  equation  (29)  for  propeller  diameters  of  2.0,  3.0, 

4.0,  and  5.0  meters.  The  curve  shows  that  on  the  axis  you  must  be  approximately 
one  diameter  away  from  the  propeller  before  a 6 dB  dropoff  per  doubling  of 
distance  occurs. 


Noise  Propagation  in  the  Plane  of  Rotation 

In  the  plane  of  rotation,  the  propeller  can  be  modelled  as  a convex 
arc  as  shown  in  Figure  21.  If  we  assume  that  point  sources  are  distributed 
continuously  along  the  arc,  the  total  energy  density  at  a receiving  point 
P in  the  plane  of  rotation  at  a distance,  d,  from  the  axis  of  rotation  is 
given  by 


2 WdS 

5“ 

c 4nr  c 
^1 


W dS 

2ttc  J 2 


S 
^1 


(31) 


Where:  W is  the  sound  power  per  unit  length  of  arc 

c is  the  velocity  of  sound 

Using  the  law  of  cosines 


r^  = d^  + R^  - 2Rd  cos  0 


r _ W f^2  Rd0 


c d +R  -2Rd  COS0 
^1 


^ J 


d0 

<5—5 — 5—^ 


^1  cos  0 


Using  the  identity  . d0 

■'a+bcos0 


(32) 

(33) 

(34) 


2 2 
0/2  ® ^ b t 

/aZ-b2 
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FIGURE  21:  GEOMETRY  FOR  A PROPELLER  AS  VIEWED  FROM 

A MEASUREMENT  POINT,  P.  IN  THE  PLANE  OF 
ROTATION 
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we  obtain 


(35) 


E = (^^<*’>)  --L 


ttC 


(36) 


?7  [-==^tan-l  (%§^)] 


TTC 


viv 


E = 


^ r 1 

TTC  '■ 


(^^)^-4d^ 


.tan  ( 


i 


2 2 

d +R  t2  ,.2 


-1 


d2.R2 


-)] 


- 2d 


(37) 


(38) 


E = 


2W  r_  1 

TTC  *■  *" 


JEE 


i 


d^-i-R^)^-4d^R^ 
tan"^  ? 5 )] 


m!e!  (d^^R^)  - 2dR 


(39) 


r _ 2W  r R ^d+R\T 

« <3Er)] 


(40) 


_ 2W 


^ r— j—  tan-l  /d/R+K-. 

wrR  ^ ? win  (h/d  I'J 

(d/R)^-l 


(41) 


Repeating  the  approach  used  to  derive  equations  (11)  through  (15),  we  obtain 

SEE  . PWL^  0-10  lo,,„*  * 10  log,„  tan''  (^)]  (42) 

Where:  PWL^^  is  the  power  level  per  unit  length 
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If  we  assume  that  the  length  of  radiating  arc  is  equal  to  the  propeller 
circumference,  we  obtain  the  following  expression  for  total  power  level 


PWL  = PWLj^  + 10  log^g2TrR 


(43) 


Therefore, 


SPL-PWL 


-20  log^QTrR  + 10  log^Q 


[ ^—5 — tan"^ 

(d/R)  -1 


(44) 


Figure  22  is  a plot  of  the  third  expression  of  equation  (44)  as  a function 
of  d/R.  The  curve  shows  that  the  plane  of  rotation  measurement  points  must 
be  approximately  3 diameters  from  the  propeller  blade  tip  before  a true  6 dB 
per  doulbing  of  distance  is  obtained.  Figure  23  is  a plot  of  equation  (44) 
as  a function  of  measurement  distance,  d,  and  propeller  radius,  R,  in  meters. 

This  result  can  be  directly  compared  with  the  result  from  equation  (6)  as  plotted 
in  Figure  11  to  determine  the  effects  of  reverberation  and  source  geometry  on 
the  noise  propagation  characteristics. 
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DISTANCE  FROM  THE  ACOUSTIC  CENTER  OF  PROPELLER  IN  PLANE  OF 
ROTATION,  d (meters) 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  has  described  the  development  of  an  indoor  acoustic  test 
facility  for  static  testing  of  fans  and  propellers.  The  performance  measure- 
ment system  has  been  used  for  many  years  with  demonstrated  accuracy  and 
dependability.  The  acoustic  characteristics  of  the  facility  have  been 
thoroughly  studied  and  documented.  The  problem  of  data  contamination  due  to 
acoustic  reflections  has  been  minimized  by  the  use  of  high  absorption  acoustic 
treatment  on  the  near  walls  and  the  floor  of  the  test  rig.  Room  mode  inter- 
actions are  not  considered  significant  above  60  Hertz  due  to  the  large  size 
of  the  room.  Reverberation  characteristics  for  the  facility  have  been 
thoroughly  evaluated  and  can  be  used  to  correct  measured  data  where  necessary. 
Finally,  background  noise  has  been  assessed  and  is  not  considered  a problem 
for  most  propeller  and  fan  configurations.  However,  background  noise  measure- 
ments are  made  an  integral  part  of  each  test  program  to  ensure  that  the 
test  data  is  not  contaminated  in  the  frequency  range  of  interest. 

In  Section  six  of  this  report,  a brief  study  was  made  to  assess  acoustic 
propagation  characteristics  with  distance  in  the  plane  of  rotation  and  along 
the  propeller  axis.  Characteristic  curves  for  on-axis  propagation  (Figure  20) 
and  plane  of  rotation  propagation  (Figure  22)  should  be  developed  for  the 
propeller  under  test.  These  curves,  together  with  the  facility  reverberation 
curves  can  then  be  used  to  select  optimum  microphone  measurement  arrays. 

For  propellers  in  the  diameter  range  of  2 to  4 meters,  the  optimum  location 
for  far- field  measurements  in  this  facility  falls  between  one  and  two  diameters 
from  the  propeller  center  in  the  plane  of  rotation.  The  farthest  microphone 
point  is  limited  by  wall  effects  and  reverberation  while  the  nearest  microphone 
position  is  limited  by  the  extent  of  the  propeller  near  field.  The  propeller 
near  field  is  generally  defined  as  the  space  within  one  propeller  diameter 
from  the  blade  tips.  Measurements  in  this  area  are  generally  free  from  inter- 
ference by  reflection,  reverberation,  and  background  noise,  but  do  not  show 
the  directivity  patterns  characteristic  of  far-field  noise.  Near-field  noise 
does  not  exhibit  radial  variation  in  sound  pressure  in  accordance  with  the 
far-field  "inverse  square  law"  where  the  level  drops  6 dB  per  doubling  of 
distance  (Reference  9), 

A number  of  useful  data  presentation  programs  have  been  developed  to  aid 
in  assessing  noise/ performance  tradeoffs  for  the  articles  under  test.  These 
programs  were  described  in  Section  V and  detailed  listings  with  sample  plots 
are  shown  in  Appendix  B.  This  plotting  capability  coupled  with  the  performance 
and  acoustic  measurement  systems  provides  a very  versatile  tool  for  conducting 
static  propeller/fan  performance  and  acoustic  tests  and  for  conducting  noise/ 
performance  trade  studies. 
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APPENDIX  A 

DETAILED  SCHEMATICS  - ACOUSTIC  DATA  SYSTEM 
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I.  NEO.  INVCIITCO 


REVISIONS 


I 


I 

APPENDIX  B-1 
PROGRAM  HORSE 
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R15,STCSA,C177777,T  20,n?0.  P77T119  HCG-tEGOR  5?74^ 

FTN,R=?. 

HAP, PART. 

ATTACH,CCAUX,CCAUX,ID=X^t43->j, 

LIRRA9Y,CCAIJX. 

LGO. 

••  ENO  OF  RECO^O 

PROGRAM  HORSF  (TNOMT, output, TAPE5  = IH=HT,TAPE6  = 0UTD|JT, PLOT) 

C •••  PPOGPAM  MORSF  IS  USEn  TO  PLOT  THRUST  VS  HORSEPOWER  FOR  VARIOUS 
C OUIFT  propeller  CONFIGURATIONS.  IT  MILL  ACCEPT  JATA  WITH  5 OR  3 

C PPM  RUNS, 

r 

PRESET  PARAMETER'?  (NA«FLIST)  I 
C lORPH  = (1)  ..  ALL  AVAILABLE  OATA  USED  IN  PLOTS 

C 2 ..  I RPM  VALUES  USEO  REGARDLESS  OF  DATft  INPUT 

C MINIMUM  THRUST  (YMIN)  = 0. 

r MAXIMUM  thrust  (V^ax)  = IfcOn. 

C minIMIJH  horsepower  (X-IN)  = 0. 

C MAXIMUM  HORSEPOWEF  (XMAX)  = 2<*0 . 

C RPM  OATA  POINTS  (WRPM(l-F))  = 60 ^ , 750  . , 9 0 0 . , 1 0 53  . , 12 3 3 . 

0 

C • CHANGES  TO  THESE  PARAMETERS  SHOULD  9E  MADE  USING  THE  STANUARL’ 

C NAMElIST  format  (EXAMPLE!  STN^UT  YMINrlO. , XMIN=13 0 . i) 

r 

INTEGER  HEAD 

dimension  HEA0(4C),  RFTA(T),  X(2,T,5),  CX(7),  CY(7),  W,<P(i(5) 
namelist/ input/  WRPM, IURP“, YMIN,YMAX,XMIN,XMAX 
OATA  MRPM/600.,750. ,90 C., 1050., 1230./ 

IORPM=l  ? YMIN=0.  « YMAX=160C.  $ X(IN=0.  i XMAX=240. 

call  plot  <0  .,1.,-T) 

1 RFA0(5, INPUT) 

YSTFP=(  Y'«AX-Y.MINI  /9.  0 
XSTEP=(X<AX-XMIN)/6.0 

C *•  PFAO  IN  NUMRF^  OF  RPM  RUNS 

2 RFAD  (5,M)  NRPM 
IF(NF.PM.fo,ij9<)9tj  go  TO  1 
IE  (NRPM.EO.O)  GO  TO  7 

CALL  REAnER(MEA0,9FTA,X, Y,NkPM,0,IORPM,NSIZ,0) 

CALL  CfnteR  (HEAO,u0) 

C *•  SET  UP  AXIS  FOP  PLOTTTNG 

CALL  AXIf  (0. ,0.,21HTHRUST  (POUNDS  FORCE ) , 2 1 , fl . , 90 . , Y1 IN , YST EP ) 

CALL  AXIS  ( 0.  , 0.  , 10MHORSEPOHFP.,-n  ,6.5, 0.  , XMIN,  XSTFP) 

ENroOF(MO, 100, HEAD)  HEAD 
lOP  FORMAT(4nAl) 

CALI  SYM90L  ( 0 . 25 , 8 . se , q . 1 5 , HEAD , 3 . , 40 ) 

C ••  SET  UP  OATA  ARRAY  FOP  CONSTANT  BLAOE  ANGLES 
90  4 1*1,3 
DO  T J*1,NSIZ 
cx( j) *y (1 ,i,j) 

T CY(J)=X(?,I,J) 

CX(NSIZfl)*XMIN 
CX(NSIZ»2) sXSTEP 
CY(NSIZ^l) *YMIN 
CY(NS1Z»?) *YSTEO 

C ••  PLOT  POINTS  And  DRAW  A SMOOTH  CURVE 
CALL  FLINE  (CX,CY,-NfIZ,l,l,  (I-D) 

best  available  copy 
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c ••  OPAW  UFGE:^0  PY  END  OF  LINE 

CALL  WLINE  (CX,CY, CNSIZ*2I  , 0.  ,XL,YL) 

CALL  SY><-»OL  ( JXLf.l)  ,YL  ,.0e,2?H—  DEG.  BLADE  ANGLE, 0.  ,22) 

4 CALL  NU-iDFR  ( { XL  ♦ . 34)  , YL  , 0.  0 0 , BET  A ( I ) , 0 . , - 1 ) 

C ••  SET  UP  DATA  AS^AY  FQP  CONSTANT  RPM*S 

00  E J=1,NSI7 
00  E 1=1,3 
CX(I)=X(1,I,J) 

5 CY( I) =X (2,1,  J) 

cx(u)=vntn 

CX (5I=XST  pp 

DY(4)=YNIN 

CY(5)=YSTPP 

C ••  PLOT  POINTS  and  hraw  A SMOOTH  CURVE 
CALL  FLIME  (CX,PY,-3,1,1, (J-1) ) 

C OPAH  LEGEND  BY  END  OF  LINE 

CALL  WLINE  (CX,CY,5,P0. ,XL,YL) 

CALL  SYMBOL  ( XL , ( YL ♦ . 1) , , 09 , 11H--  RPM,O0.,ll) 

C ••  COMPUTE  DESIRED  INDEX  FOR  TITLE'S  RPH  VALUE 
LABFL=IFIx(FLOAT  (J)  *(e:.0/FLOAT(NSI2)  M-.OOl) 

9 CALL  NUMBER  (XL, (YL».34) ,.CB,WRPM{LA3EL) ,90. ,-l) 

CALL  PLOT  (9.«;,f|.,-j) 

GO  TO  2 
7 CONTIN'JP 
call  PLO'^E 
STOP 

j 

9 FORMAT  (1015) 

FND 
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SUBR0UTI*^E:  ?FanpP  (HrAO.PFTa^X.Yj^JRPMjNHIKEflOf^PH.NSTZjTOjjC) 
3INENSI0M  HEA0(4"),  X(2,3,5),  Y(U,3,';), 

INTFGfr  head 
REAT  (5,14)  HFAn 
IF  (lOPPH.FQ.l)  r,0  TO  9 
IF  (I0RRM.En.?.AS3,NPPK.El.?l  GO  TO  9 
C READ  3 POINTS  TF  5 APE  AVAILABLE 
30  P 1 = 1,  3 
REAP  (5,15)  5FTA(I) 

00  1 J=l,2 

PEAO  (5,16)  »"M,X(1,I,J) ,X(?,T, J) 

1  REAP  (f:,!?)  TRA«;h 
? READ  (5,16)  ko*<,X(l,T,3)  ,x  ( ’.T,  3) 

8 NSI7=x 
GO  TO  13 

C ••  READ  ALL  AVAILAaL*"  DATA 

0 on  10  1=1,3 

READ  (*:,15)  PFTA(I) 

00  10  J=1,NRPM 

10  READ  (5,16)  RP‘<,X(1,I,J)  ,X(2,I,  J) 

MSI7=N?»-1 

13  RETURN 

0 

14  FORMAT  (40A1) 

15  FORMAT  (T0X,Fin.3) 

16  format  (3F10.T) 

17  format  (P10.3) 

ENO 

f 

f 

SUOPOUTINE  CFNTpp  ( HFAD,L'=‘NGTH) 

C SUBP0UTI>(E  OENTer  is  U560  TO  CENTER  A HEAOTNG  WITHIN  A VA’lAfL'^- 

C length  title  LINr 

INTEGER  HFAO 
DIMENSION  Hf«0(LFNGTH) 

IM<  = 0 

C *•  COUNT  THE  «IGHT-maN0  ^LANKE 

1 LOC=LENGTH-In< 

IF  (heao(loo)  .nf.ih  ) ro  '•o  ? 

TNK=IN<»t 
GO  TO  1 

C ••  PLACE  HAL=’  THE  OPTECTEO  BLANKS  ON  THE  HEADING’S  L"P’T  iID“ 

2 MOVFsIFIx (( (FLOAT (INK)) /2) ♦.rCl) 

IF(MOVE.FQ.n)  RETURN 
00  4 I=1,M0VR 
N»LFNGTH-1 
IP0E=LFNGTh*1 
00  ■*  J=1,N 

3 HEAD(IP0S-J)=HEA0(LENGTH-J) 

4 HEA0<I)=1H 

return 

ENO 


%“>  o 


SUfJPO'JTI.'JF  WLINF  (X,V,N,  W''H,  XL,  VL) 

•••  SUBPOUTI'JE  WLINP  flPRflNGES  E'aTfl  SETS  IN  DECREASING  O'^OEK  ^ESEECTIVE 
TO  FITHr.-*  X OP.  V.  IT  THEN  PETURNS  THE  PHYSICAL  PAGE  CODROINATES 
or  THE  DATA  SE^  FOUNP  GREATEST  RELA^tve  TO  THE  OTHERS. 
dimension  X (N) , Y(N) 

IP?=N-3 
lP3rN-P 
DO  T 1=1, IP3 
IPl=I»l 

00  3 J=IP1,IPT 
IF  (MCH.ro. DO.)  GO  TO  1 

C ••  COMPAPr  "X"  VALUF> 

IF  (X(I)  .GE.X(J) ) GO  TO  3 
GO  TO  ? 

C ••  COMPAPE  "Y"  VALUES 

1 IF  (Y(I) .GF.Y( J) ) GO  TO  3 
C ••  interchange  the  OUT-OF-OROrP  X-Y  OATA  SETS 

? TpMprX(I) 

X(I)=X( J) 

X( J)=TENP 
TENF=Y( I) 

Y(IJrV( J) 

Y( J)=TEM? 

3 CONTINUE 

C ••  COMPUTE  THE  PHYSICAL  PAGE  COOROINATES 
YL=  (Y(l)-Y(N-D)  /Y(N) 

XL=  <X(1 )-X(N-l>) /X(N) 

RETURN 

FNO 

- END  OF  PEPORO 


US!  ®'' 
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I 


5INPUT 


•5 


f- 

1 16.75 

9.33 

60?. 

17. 

58. 

7*;o, 

??. 

93. 

ROl. 

30. 

134. 

1050  . 

33. 

186. 

1?00. 

5L. 

?4«;. 

6 

1 16.75 

9.33 

601. 

?6. 

138. 

^5?, 

37. 

?96. 

qoo. 

56. 

430. 

1050  . 

85. 

601. 

1?01. 

1?1. 

736. 

6 

116.75 

9.33 

600. 

36. 

34?. 

751. 

61. 

537. 

900  . 

104. 

785. 

1051. 

162. 

1081. 

i?ca. 

240  . 

1430. 

99999 

SINPUT  X-flXs 

120.,  Y**ax  = 30D.* 

3 

? SETS  (6"SHRT,3"aiHRT, 

NORM) 

6 

630. 

8. 

?. 

90?, 

21. 

5. 

1?00. 

39. 

11. 

599. 

9. 

62. 

899. 

25, 

117. 

1?00  . 

49. 

?17. 

6 

599. 

14. 

155. 

899. 

<•5. 

3se. 

1?00  . 

08, 

658, 

— 9LAN< 

Qlie't 

5.0  4,0  1.0153 

10.1  4,0  1.0133 

15.0  4.0  1.0133 

1.1 

5.1 

10.1 
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6 BLROES.  STRNOflRD  o 2 SETS  ( 6"SHRT. 3“SHRT. NORM ) 


Figur«B1-1:  Plots  From  Sample  HORSE  Input  Data  Cords 


APPENDIX  B-2 


PROGRAM  THORSE 


^•16fSTCSA,C‘177'’77,T4P,IO<<0.  P72011<j  MCGREGOR  55421 
PTN,R=3. 

MAP, PART. 

ATTACH,PL0TLI5,  OGIVFFILE’,CY  = 2,I3=’7  40  03!3  ,SN=AS0. 

ATTACH,CCAUX, CCAUX, ID=X654321. 

LIBRARY, CCAUX. 

L0A0,L50. 

SATISFY. 

L3SET,LIB=PLOTLI9. 

SATISFY, PLOTLIT, 

EXECUTE. 

••  FMO  OF  FECORO 

PROGRAM  r HORSE  (INPUT,  OUTPUT,  TAPES  = INPUT,  TAP E6  = J'JTPUT, PLOT) 

C •••  THORSE  IS  USEO  TO  PLOT  T riRUST/HORSEPOWER  VS.  HORSEPOWER  FOR 
C VARIOUS  TUIET  PRQPELLOR  CONFIGURATIONS.  IT  WILL  ACCEPT  OaTA 

C with  3 0?  5 \PM  POINTS.  ALL  OR  PART  OF  THE  DATA  MAY  BF  USED. 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


PRESET  PARAMETERS  (NAMELIST)  t 

lORPM  = (1)  ..  ALL  AVAILABLE  DATA  USED  IN  PLOTS 

2 ..  3 RPH  VALUES  USED  REGAROLFSS  OF  DATA  INPUT 

YMIN  = (3.)  ..  MINIMUM  T HPUST/HORSEPOWER  FOR  PLOT 

YMAX  = (11.)  ..  >-AX1mUM  THRUST/HORSEPOWER  FOR  PLOT 

XMIN  = (0.)  ..  Mir.'IMUM  HORSEPOWER  FOR  »LOT 

XMAX  = (2*-n.)  ..  MAXIMUM  HORSEPOWER  FOR  PLOT 

RPM  DATA  POINTS  (WRFM(l-E))  = S Q 0 . , 75 0 . , B 00 . , 1 0 pO . , 12 J 0 . 


C * CHANGES  TO  ’’HESF  PARAMETERS  SHOULD  BE  MADE  USING  THE  STANuARD 

C NAMELIST  FORMAT.  (EXAMPLE)  TINPUT  YmI N=6 . , YM A X= 24 . S ) 

C 


INTEGER  HEAD,OFD 

DIMENSION  H£AD(40),  °ftA(3),  X(2,3,5),  CX(7),  CY(7),  WROm(5) 
1,XC(102),  YC(102),  ORD(o),  COEF(^) 

NAMELIST/ INPUT/  WPPM,IORPM, YMIN, YMAX, XMIN, XM AX 
OAT  A WR?«/S00  ., 7^0., 900.,  lO'^O.,  120  0./ 

IORPM=l  ' YMIN=3.  S YMAX=ll.  $ XMIN=0.  ? XMAX=240. 

CALL  PLO^  (0.,1.,*3) 

1 REAu(5, INPUT) 

YSTEP=(  YMAX-YMIN)/#.!^ 

XST«-P=(X1AX-XMIN)/6.P 
C **  READ  IN  NJ“BER  OF  RPM  RUNS 

2 READ  (5,5)  NRPM 

IF  (NRPM. to. n)  GO  TO  7 
IF(NRPM.EO.B9BBB)  GO  TO  1 

CALL  REAjtR  ( HE An, q£T A , X , Y , NRPM , 0 , lORPM, NSIZ , 0 ) 

CALL  CENTER  (HFAC,40) 

C ••  SFT  UP  AXIS  FOR  PLOTTING 

CALL  AXIS  (0 .,0. ,17HTMFJST/HORSEPOWER,17, 5. ,B0. ,YMIN,YSTEP) 
CALL  AXIS  (0. ,0.,10HWORSE“OWER,-13,6.5,0.,XMIN,XSTEP) 
ENCODE(40, 100, HEAD)  HEAD 
100  FORMaT(40A1) 

call  symbol  ( 0.50,'*.  35,  a. 15, head, 0.  ,40) 

C *•  SET  UP  DA'^'A  AK?AV  FOR  CONSTANT  BLADE  ANGLES 


00  4 Isl,  3 
DO  3 J=l,NSI7 
CX( J):X (1,I,J> 
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3 CV (J) =X  (2, 1 , J) /Cv  ( j) 

if(cy (Njsiz)  ,LT ,Y /IN)  ro  u 

: ••  O-^AW  A SN'IDTH  CU!VE  A^'j/QC  T^TA  POIMT3 

CALL  CJ-^\/FT(CX,CY,NSIZ,XMT,4,XST£",Y.iIN,YCTF3,-?,l,(i-l),l.,0.) 
C ••  DRAW  LEGn/3  3'^  t OF  LlN«^ 

XL=(C< (WCiZ) -XNIN) /XFTtP 
VL=  (CY  ( )£IZ)-Y‘1In)/Y'^1FP 

CALL  SY^’OL  ( ( XL+- . 1)  , YL,  .05  , ?2H—  3 EG  . 6LA1E  AwGL'1 , 0 . , 22) 

CALL  NJN3EP  ( (XL«-.34)  ,YL,a.T6,  JE^Ad)  ,o.,-l) 

4 CONTINUE 

C SET  UP  3A*A  ARRAY  Fnp  CONSTANT  RP’I'C 
00  o J=1,NSIZ 

00  ‘p  1 = 1,1 

CX( I)=X (1 , I,J) 

5 CY ( I) =X {’ , I, J) /CX (I ) 

Ms  Y 

IF(rY(l) ,GE. YMIn)  GO  TO  50 
00  »n  <=1,4 
CX(<)=Cx(t^*l) 

40  CY(<)  =CY(K4-1) 

cx(  .3)  = n. 

CY ( 3) =n  . 

M=2 

C •*  PLOT  POINTS  ANO  OkAW  A F>103TH  CURVE 

5 0 CALL  CU?VFT  ( C X , C Y , M , XMIN  , XST  F P,  YMI  ij , YS  TE  P , -2 , 1 , ( J-1)  , 1.  , 0.) 

C DRAW  LEGEND  OY  EN3  OF  LINE 
XL=(rX('1)  -XMI  >1)  /XCTEP 
YL=  (CY  ( i)  -YMIN)  /YSTFF 

call  symbol  (XL,  (YL»-.1)  , .3S,11H—  Rpm,:»0.,11) 

C ••  COMPUTE  JESIRtO  IMOL'y  FOR  TITLE'S  R^M  VALUE 
LA  pels  I FIX (FLO A* ( J)* ( E . 0 /FL 0 A T ( NS I Z ) ) +.3  01) 

6 CALL  NUNJER  ( XL  , ( YL  <• . 3 4)  , . 0 3 , WRPM  ( L A-3EL)  , 00  . 1 ) 

call  pLO'  (9. 5,0., -3) 

GO  TO  2 

7 CONTINUE 
CALL  PLO'L 
STOP 

C 

3 FORMAT  (1015) 

ENO 
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C'  o 


SU'^-'OUTI'Jt  {HEAD,  3«- T a , X , Y , Nill  KF  , IO,^p  >1,  NSI  Z , I 0 5 ) 

3IMFN5I0N  HFfi[)(‘*C),  X(2,3,5>,  ‘?FTA(3> 

INTEGER  HEAD 
pEAO  (E,i4)  Mean 
IF  ( lORPM.EO.l)  GO  TO  9 
IF  (IORP-.EQ.Z.ANO.NRdm.FQ.T)  go  to  9 
C ••  KEAO  3 POINTS  IF  9 ARE  AVAILARLE 
00  2 1=1,3 
RFAO  (5,15)  BF-^A(I) 

00  1 J=l,2 

READ  (5,16)  RPM,X(1,T,J) ,X(?,I, J) 

1 REAJ  (5,16)  TRA5H 

2 REAO  (5,15)  PPM,x(l,T,3) ,X(2,I,3) 

*1SIZ  = 3 

GO  TO  13 

C ••  READ  all  AVAILABLE  DATA 
9 00  in  1=1,3 

REAP  (5,15)  B'^TA(I) 

DO  10  J=1,NRPM 

10  READ  (5,15)  PPH, X(1 , I, J) ,X (?,I, J) 

NSI2  = NRP-3 

13  RETURN 
C 

14  format  (40A1) 

15  FORMAT  (30X,F10.3) 

16  FORMAT  C'FIO.3) 

FNO 


SU9R0UTIMF  CENTER  ( HEAD , LENGTH) 

subroutine  CENTER  IS  USEU  TO  CENTER  A HEADING  WITHIN  A VARIABLE- 
LENGTH  TITLE  LIN«^ 

integer  head 

OIMENSION  H£B  0 (LENG'fH) 

INK  = C 

C ••  COUNT  THE  right-hand  BLANKS 

1 L0C=LENGTH-INK 

IF  (HEAO(LOC) .NE.IH  ) GO  TO  2 

INK=INK+1 

GO  TO  1 

C ••  PLACE  HAL-  THE  DETECTED  9LANKS  ON  THE  HEADING'S  LEFT  SIDE 

2 MOVF=IFIY ( ( (FLOAT (INK)) /2) ♦.001) 

IF(HOVE.EQ.O)  RETURN 

DO  4 1=1, HOVE 
N=LENGTH-1 
IP0S=LENGTH*1 
DO  3 J=1,N 

3 HEAn(lP03-J)=H£AO(LENGTH-J) 

4 HEAD(I)=1H 


RETURN 

FNO 

END  OF  FECORO 


M---'  COVX 

h»‘-'  — 
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3 


SIMPUT  S 

5 

0 3LA0ES,  STAN)ART 

116.75 

9.33 

5.3 

4.  0 

1. 0133 

S02. 

13. 

58. 

750. 

22. 

93. 

901. 

3C  . 

134. 

1050. 

36  • 

196. 

1200, 

54  . 

2US. 

f 

116.75 

9.33 

10.0 

4.  0 

1.0133 

501. 

26. 

198. 

752. 

37  . 

298. 

900  . 

56. 

430. 

1050. 

9-^  . 

e 01. 

1201. 

121. 

796  . 

6 

116.75 

9.33 

15.0 

4.  J 

1 . 01 35 

oOO. 

36. 

34?. 

751  . 

51. 

537. 

900. 

1 0^. 

78p. 

1051. 

162. 

1091. 

1200  . 

2**C  . 

14  3 0. 

99999 

?INPUT  XMAX= 

120. » 

3 

2 SETS  (6“?HRT, 3"SHRT, 

NORM  ) 

6 

1.0 

63  0. 

5. 

2. 

902. 

21. 

5 . 

1200  . 

39. 

11. 

6 

5.3 

599. 

9. 

52. 

899. 

25. 

117. 

1200. 

h9  . 

217. 

10.3 

599. 

lu. 

155 . 

899. 

45  . 

356. 

1200  . 

9e. 

556. 

BLANK 

CARO 

Ril,STCSA,Ci5i3  00  ,T20,I02»).  P''2niy  MCG  >EGOP  527W-+ 
FTN,R=3. 

MAP,PA«?T. 

ATTACH,CCAUX,Ck.AUX,  ID  = Xbb*.321. 

LIBRARY, CCAUX. 

LGO. 


END  OF  RE:0^0 


“ROGRan  ►JOlSYl  (INPUT,  OUTPJT  ,TAPE5  = InPJT  , T AP£r3  = oUTr  UT  , PLOT  j 
C 

C PROGRAM  MOISYl  WILL  PLOT  THt  SOUUO  o<ES3'JR£  LEVEL  VS.  F^_(jU«lNv.Y 

C PLUS  THE  OVERALL  FOR  ANY  PROPELLOR  RUNS.  THE  U^E-‘.  MAY  SP-JJIFY 

C THE  OATA  CARD  FORMAT,  THE  OUTPUT  PLOT  RANGE,  AND  WHi:HF.>  J<  NOT 

C HE  DESIRES  A DBA  CURVE  FOR  EACH  SPECTRA  plOTTEJ.  DEFAULT  V,.LUcS 

C ARE  LISTED  BELOW  WITH  INST^UoTIOJi  FOR  THEIR  CHA.JG'".  GAIA  W.lc 

C 3E  varying  PPM  RUNS  FPC-H  A SINGLE  ITCRUPHONE.  (FOR  A SING*.E 

C RPM  FROM  SEVERAL  MICROPHONES,  USt  N0ISY2) 

C 


C PRESET 

r 

PARAMETERS  (NAMELIST) 

r 

otform 

= (0) 

3.15  THROUGH  20K  PLUS  OVERALL 

c 

r 

1 .. 

20 

THROUGH  2Q< 

PLUS 

UVERALu 

o 

c 

I03A3T 

- (0)  .. 

NOM'ai.  PLOTS, 

NO  ADPITIOrVA^  DBA 

CU-IVE 

c 

p 

1 • • 

AOOniONAL  U5A 

CURVE 

DESIRE J 

c 

SPlMIN 

= (20.) 

• • 

MINIMUM  SPL 

VALUE 

FOr  Pi-OT 

c 

p 

SPlMAX 

= 1120.) 

• • 

MAXIMUM  SPl 

value 

Fur  HwDT 

c 

start 

= (1)  .. 

starting  frequency  INDEX 

c 

p 

finish 

= (2tJ).. 

CLD 

SING  FFEUUENCY  INDFX 

c 

INOEX  TMPLE 

c 

I HZ 

I HZ 

1 

HZ  I 

H7 

I HZ 

I HZ 

c 

1 20 

b 63 

1 1 

20  0 16 

63  0 

21  200J 

2o  j30: 

c 

2 25 

7 <*0 

12 

230  17 

90C 

22  2533 

2 7 i 0 0 D 

c 

3 32 

6 ICO 

1 3 

313  IS 

lOOG 

23  3150 

23  IQOCO 

c 

4 43 

0 12E 

14 

hOO  19 

12=0 

24  4UJ3 

c 

? 50 

10  15C 

15 

50g  20 

IbtiC 

25  50G0 

c 

C • CHANGES  TO  THESE  PARAMETERS  SHOULD  HE  HADE  USING  THE  jIANDL-.L 

C namelist  FDR-AT  (EXAMFLL*  iINP'UT  FINISH=13,  IuOA5T  = ls) 

C 


INTEGER  DTFDRM,START,F1NI3H 
OIHENSION  XFR (31) 

DIMENSION  TITLEX(h)  ,TITL£Y  ( - ) ,A  ( Jl>  ,FR{  30)  , TITLc(h)  ,UNX(  -»>  ,UNY  (-,) 
NAMELIST/ INPUT/  uT  FO  A rl,  I 09A  ST  , ST  ART  , PI  NI SH  , SPL  M AX  , SPLMIN 
DATA  TITlEX/IOHONP  THIRD  ,IOHOCTAVE  CEN,10hTE?  FhEQUE, 
llOHNCY  IM  HZ  /,TITLEY/10H50UMD  LEVE,10HLi  IN  DH  R,10HL  0.0002  M, 
210HICR0SAR  /,Fp/2  0.  ,25.,  31.S>,40.,SO.,b3.,30.0,100.0,l2b.a,lDr. , 
320  0. * 2b0., 31b. , 40n., 500., b 30. , dOO.,  1000.,  12i>0.,lbQG.,20CC-., EEOC., 
•»315  0.,H00  0.,e0U0.,b33  0.  , 90  00., 10000., 10., 3. 4<U/ 
5,UNX/4.0,&.2b,0.,l./,UNY/1.75,1.75,0., 1./ 

OTFOPMsO  I IC'3ASTsO  i START  = 1 S FlNlGH  = 2d  S SPLMAx=120, 

SPLMIN»20. 

CALL  Plot  (0.,  1., -3) 
i REA0(5, INPUT) 


BEST  AVAILABLE  COPY 
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1 


1 


r 

n A J ( 3 , l : C ) 1 I T L r , I c 

1^  n?"!  FOR’^31  (-ttiO,in) 

r IFdTfi  nC  1)  .LQ.  lOHNtWPAiRflM  ) GO  TO  5 

I TF(TITLE(  1)  .EO.lOHPNrUATA  >GO  '^O  <^99 

' HM=STftRT 

NN=FINISH 

I CHNn=ICHMO 

SPLSTFP=(  SPL»1flX-9PLHTH» /8. 0 

11=0 

YLOCsl.SS 
NUMPTSsNN-M'I^T 
'^••OOAM  CNO  LA9FL  btfS 

CALL  AX  19  ( 0.  ,0.  ,TI'^LFY,4n,e.,90.,SOLMlN,SPLSTEP) 
tjnCs&.G/FLOATJNU'IPTS-I) 

TIC<=0 

CALL  PLOT (6.0 ,p.n,X) 

CALL  '’LOT  (0.0,  0.  0,2) 

N=1 

DO  too  IOT=MM,NN 
IFdPT.GT  .2)N  = -1 
reuL  PLOT (TICK, -0.1, 2» 

CALL  NU'' JFR(Tir<,-0.  1,  0,  Qt  , FF  ( I PT  ) , -60 . 0 , N) 

TICK=TIC<»TIA'C 
100  CALL  PLOT  (TICK, 0.0, 31 
CALL  PLOT (6. 0,0. 0,3) 

CALL  P( OT(6.0,-3.1,?) 

I CALL  SY1POL(b.35,-0.2,0.1,3HO/A,0.0,-») 

CALL  SYMIOL  { .75,-.6,  .1,44H0N£-THIRD  OCTAVE  9A^^D  CENTER  FREQUENCY  I 
2N  HZ,0.0,44) 

CALL  9 YHPOL  (1.35,6.  , 0 . 1,  T TTLE , 0 . , 40) 

CALL  SYM3oL(?.0F,8.5,0.15,11HVIKE  NUN9ER , 0 . , 11 ) 

CALI  NU'^OER  ( 3. '»5,6 .5,0.  15,  CHNO,  0.  , -1) 

CALL  SY*10nL  (4.5,1.'^tc;,o,i,  3hrPM,0.,  T) 

CALL  LINP (UNX,UNY,?, 1,0,  0) 

C**PUT  DATA  IN  A ARRAY 
200  RpaD(5,  l^ODIRUNNO 
RUNNO=IP'JNNO 
ID9A=in9‘'ST 
1001  FPRMAT(IIO) 

IFdPHNNO.EO.niGO  TO  899 

CALL  GETTN(DTF0RK,1RUNN0,ICHN0,A,H'1,NN) 

C»*PUT  IN  SYM90L  PQR  OVERALL  SPL 
( 1002  OALL=OVE?AL  (A,M-1,NN) 

WRITF(6,2000) OALL 
2000  FORMAT (PIO. 4) 

OALLs(OALL-SPLMIN) /5PL5TEP 

CALL  SYHPQL (6.0 0,OALL,C.l, 11,0. 0,-1) 

CALL  SYM10L(4.35,YLOCf.05,0.1,II,0. ,-l) 

CALL  NUf1TFR(H.7,YLOf ,0.1,PUNNO,0.,-1) 

IF( iniAST.EO. O.OR.IDCA.EO.l)  GO  TO  249 
CALL  SV9J0L (E ,2S,YLOC,0. 1,5H(DBA) ,0. ,5) 

c** remove  out-df  range -points 

249  DO  3f0  JJ=1,KUMPTS 
IF  (A (J J) .GE.SCLHIN)  GO  TO  250 
A( JJ)=?.l25fSPLlIN 

250  IF  (A( JJ) .LE.SPtNAX)  GO  TO  300 
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I 


4( JJ)=5P.MAX-i.l?5 
100  ^ONTI'JJE 
C**0RAW  EACH  LINE 

A(NUHPT5»1) =SPLhIN 
A (NUMPTS»?)  =iPL"^  i EP 
XFR(NJHPT  S*l) =0 
XFRCNUHPTS^Z) =1. 

JJ  = HH 

DO  ?0l  I=1,NUMPT«: 

XFR(I) =TINC*FL0AT(I-1) 

201  DONTINJE 

CALL  LIN-:  (XFP,A,NU“PTS,l,lf  ID 
YLOC  = YLO:-O.l'5 

IF  (I33A.HE.1)  GO  TO  799 
CALL  03Al1N(A ,mh,nN) 
I09A=ID3A-1 
GO  TO  1002 
799  CONTINJE 
GO  TO  20? 

899  call  PLOT(8.5,O.C,-T) 

GO  TO  10 
999  CONTIN'JE 
CALL  PlOTE 
STOP 
END 


SUBROUTINE  GtTIN (<, IRUNNO, ICHNO,A , MM,NN) 
DIMENSION  A(31) 

IF  (K.EQ. 1)  GO  To  50 
REA0(9» 1000) ( A(L) ,L=1,2) 
READ(5,1001)(A(L),L=3,12) 
REAO<5flOOl)(A(L),L=13,22) 

READ(3>,  1003)  (ACL)  ,L  = ?3,2S) 

1000  F0RMAT(62X,2F6.1) 

1001  FORMAK  13F6.1) 

1003  F0RMAT(SF6.1) 

GO  TO  75 

50  READ  (5,2000)  (A(L),L  = l,n) 

READ  (5,1001)  (ACL) ,L:ll,20> 

READ  (5,2001)  (A (L) ,L=21,26) 

REAP  (5,2002)  GARRAGr 

2000  format  (16X,10F5.1) 

2001  FORMAT  (iFS.l) 

2002  F0RMAT(F3.1) 

75  MRITE  (5,200) 

MRITE(6,10  0)  (A(l)  ,L>Hrl,NN) 

100  FORNATdC  ,10E12.2) 

200  FORMAT  (2X, 10HO9 


RETURN 

FNO 


BBI.AVMIASLE  COPY 
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SUS^OOTUE  OjftLlN(A 
DIMENSION  A (311  , D'^AWT  (23) 

JATfl  J94U/-tO.'?, 

A“44»7f*3)*4>*34»&f“30*2f'"2o*2>*22*5>*19*l>*16«lf“lJ**+>*10«Jf“d*o> 
B“6»6>“4*3>**3«2f*l»9f“n*3f0»nj0»6>l*0fl*2fl«3>l«2>l»0f3*5>“[)#lf 
C-1. l,-2.>/ 

30  100  1=1,26 
100  A(I)=A(I)  ♦03AWT  (I) 

WRITE  (6,  200) 

WRITE  (o,  300)  (A(L)  ,L  = HH,»4N) 

200  FORMAT  (?X,10H03A 
300  FORMAT  (IH  ,10E12.2) 

RETURM 

ENO 

••  END  OF  Ri:o^D 


~U-4CTIOM  OOErsAL  (A,m-,nn) 
JIMENSIO-J  A(J1) 

SUM=0.0 

DO  100  I=“‘1,NS 

100  SUMsSOMK  10.C*»(A(I)/10.0)  ) 
Ol/Er  Al  = 1J.0*ALOG10(SUM) 
RETURN 
END 


fINPUT  OTFORl=l,  I09AST=1* 

6 BLADE  BASELINE  - IE  CEGR'^E  PITCH  3 

ROO 

:DI  X*R003  C003*  060.7  054.5  037.5  062.2  064.7  063.2  07^.2  082.5  066.2  073.0 
075.7  075.5  0^4.2  079.5  079.5  Obl.O  030. » 080.5  091.6  090.7 
081.2  08  1.5  ObL.O  081.  0 081.7  03^.0  03>i.  7 083.  0 062.5  048.7 
044.0  094.5 

— blank  CARO  - — 

ENOOATA 


BtSl,  AVAILABLE  COPY 


70 


ONC-THIRO  OCTRVe  WWO  CCNTCR  FREQUENCY  IN  HZ 
Figur*  B3-1:  Plot  From  Sompio  NOISYl  Input  Data  Cards 


APPENDIX  B-4 
PROGRAM  NOISY2 


•^AC,STCSA,CMo0'J00,T'»0,  lose.  P72T119  mCG It r,r)K 
PTM. 

ATTaCH,CCKUX,OCAUX,  IC'  = Xo^  US21. 

Ll3'?AH;y  ,ccajx . 

LGO  • 

••  END  OP  Kr:j-;D 

PkO&r<A^  '101'=;Y2  (Input  ,0HTPJI  ,T/'Pt&  = I.J  = JT  , T SPEo  = OUT  Pu  f ,PLO  i ) 

c 

C PROGRAM  40I<5Y2  WILL  •’LOT  SQ'JNJ  P'ESSJ'^E  LE«/Zu  \h, . P.<^  Jj[  f<'.:  Y 

C »LUS  HE  3«/£R/4LL  FOP  A^gy  ^kOFElLO?  t.UNb.  ’H£  UiE,<  MAY  -Hc.ZIFY 

C the  TATA  CARD  FO-Mai,  THfc  OUTC’IT  -mUGL,  AfJJ  WHLHE<  JR  f.ji 

C HE  OESKES  A UOA  CURVE  POr  -ACH  SPPvTRA  PLOTT'j.  uFFAOlT  VAlJcJ 

G ARE  LiSTtJ  3ELOW  WITH  INSI  RUCTIONS  POP  Htl<  CH.ggGi..  DAT',  ki.t 

C IE  A :ON:>TaNr  PPM  PROM  ANY  igUHOET.  UP  MIGfs  JPWO^cS.  (FOR  fl 

C single  MiCROPHOWP  AMP  VARYInG  RP'I'S,  USE  NOISYl) 

r 

Z ****  Pf'ESt^  PA'-AHETtPS  (NflrttLiST) 

C 

C PTPOR-'  = (0)  ..  3. IE  HKOUGH  23<  =LUS  OVERALL. 

c 1 ..  20  Through  20<  pujs  overall 


I03AST  = (0) 
1 


..  NO^mAj.  PtOTSt  '10  AOUITIUUAL  LSh  CU<VC 
..  AbJlTXON'.L  03*.  CUEVF  JESIREJ 


SPlMIN  = (2n.)  ..  PINIrtUR  SPl  VALUE  FOR  PcOT 

Spl-AX  = (120.)  ..  ‘’AXIPUP  JPL  VALUE  FOR  PLOT 

STAtT  = (1)  ..  STmRTING  frequency  ijdex 
FINISH  = (26)..  CLJEIN}  FPEOUEUCY  INDEX 


INOEX  TA3LE 


c 

I 

HZ 

I 

H7 

I 

HZ 

I 

HZ 

I 

HZ 

i 

HZ 

r 

V 

1 

20 

6 

63 

11 

20C 

16 

t.30 

21 

III  DO 

2j 

j3(j  0 

1^ 

2 

25 

7 

SA 

12 

2:>0 

17 

300 

22 

250  j 

27 

«jnr, 

C 

3 

32 

8 

100 

13 

315 

18 

10  JO 

23 

il5u 

23 

1 j 0 J 0 

c 

4 

40 

9 

I’E. 

14 

■«ao 

19 

1260 

2-* 

•4  0 DQ 

c 

5 

50 

10 

150 

15 

»0U 

20 

16J0 

25 

50  JO 

• CHANGES  TO  THESE  PARa^ETEKj  SHOULD  UE  nAO<£  JSInG  THE  aTAi.JA^i 
NAME. 1ST  FOR^^AT  (EXAMPLtI  $1i(PUT  FINXS'H  = 13,  IlTiAST  = 1;) 

INTEGER  OTFukM, start , FINISH 
OIMENSION  XFP.(31) 

OIMENSIO'4  TITLEX(h)  , Til  ley  (■*),  a (31)  ,Pr(  30)  , title  (h)  ,U.<X(>)  ,UOY(n) 
NAMFLISTXINPUT/  OTFOkM,IJ3AST,START,FINISH,SPwMAX,SPlMIN 
DATA  TITlEX/IOHONE  THiRj  ,10HOCTAVE  CEN,iOHTE<  FRlQIJii, 
llOHNCY  IN  HZ  /,TITLFYY10HSOUNO  LE7E,1GHlS  IN  DO  R,  10‘(E  0.00P2  M, 
210HICRD3AR  / , FR/2 0 . , 2 5 . , 31.5tA0.,5n.,63.,30.0,100.0,12v.J,lv0., 
3200.,25  0.  , 316.  ,AOO.  I vOO.  tS31. » SCO.  f IGOa. , 12ia.  , loOO.  ,20  0 V.  ,2-iUO.  , 
4315  0., 400  0. ,5  00  0 .,6300.,  600  0., lOOaO.,10. , 3.431X 
•, UN  X/ 3. 9, 5.65,0. , 1 . / , UNY/ 3 , 0 , 8 . 0 , 0 . ,1,/ 

DTFOkMsO  ? ID1AST=0  « STA^T^1  i FINTSH=2S  * SPL1AX=123. 

SPLMIN=20. 

CALL  PLOT  (0,,  1 . ,-.3) 

5 READ  (5,  INPUT) 

10  REAr)(5,10C0)TI^LF,IRPHO  HrrT  Alt  ah  . 


BKrjVAIUSlE  COPY 


73 


1000  FORMAT  ((*A10,Iin) 

IF(TTTLE(  1)  .FQ.  lOMMFWPARAM  ) GO  'O  <5 
TF(TITLF(1) .FO.IOmENDOATA  )G0  to  990 
9PH0=IROM0 
I03A  = I03  fST 
11=0 

rLOC=7.-^5 

MMsSTART 

NN=FINIS9 

NUMPT9=MN-MM»l 
C*»ORAW  ANO  LA'iFL  AXES 

SPLSTEP=(SPLMAX-fpLMIM) 0 

CALL  AXIS  ( 0.  ,0.,TITLEY,4  0,«.  ,9''.  ,SPLMIM,S-»LSTE«») 

TINCsS.O/FLOA"^  (NUM=T«--l» 

TICX=0 

CALL  PLO^ (6.0, 0.0, S) 

CALL  PLO’-(0.0,0.n,?l 
N=1 

00  100  IPT=MM,NM 
TFdOT.G''  .2»N  = -1 
CALL  OL0T(TIC<,-0.1,?) 

CALL  NUM3ER(TTCK,>C.1,0.07,'^P(IPT)  ,-60.0,N) 

TICK=TICK+TINC 
100  CALL  PLOT (TICK, 0.0, 3) 

CALL  PLOT (6. 0,0. 0,3) 

CALL  PLO-r  (6.0, -0.1, 2) 

CALI.  SYM9  0L(6.8P,-0.2,0.1,  3H0/A,  0. 0 , ■») 

CALL  SYM90L(.79,-.6,.l,U4H0ME-THIRn  nCTAVE  OAMO  CEMTER  F^cQUErJCY  I 
?N  HZ, 0.0,44) 

CALL  SYM10L(0.0,0.5,0.15,TITLE,0.0,40  ) 

CALL  StnOL  (3.9,7.75, 0.10,  IIHMIKE  MUMOER,  0 . 0 , 1 1 ) 

CALL  SYMJOL  (4 , 875 , S . •’25,  0.  1 , 3H9RM,  O . , T) 

CALL  NUMBER (4. 37 5, 4. 026 , 0. 1 , PPMO, 0 . , -1 ) 

CALL  LINE(UNX,UMV,2, 1,0,0) 

C**PUT  DATA  IN  A AORAY 
200  REAO(5,1001)I-I<0 
XMIK0=IMTK0 

1001  FORMATdtO) 

IFdMIKO.EQ.OI  GO  TO  399 

CALL  GETIN(0TF0PM,I©UNN0,ICHN0, A,MM,MN) 

C*»PUT  IN  SVM90L  FOR  OV«=‘RALL  SPL 

1002  OALLsOVERAL(A,mm,NN) 

NRIEE(5,?000) OALL 

2000  FORMAT (Fin. 4) 

OALL*(OALL-SPL*«INI  /SPL^TE® 

CALL  SYM90L(6.00,OALL,0. 1,11,0. 0,-1) 

CALL  SYMBOL (5.5,YLOC».D5,0. 1 0,11, 0. 0,-1) 

CALL  NUM9ER(5.1,YLOC,0.10,XMIKO,n.0,-l) 
IFdOaASY.EO.O.OR.IOPA.EO.l)  GO  TO  249 
CALL  SYM9OL(9,75,YLOC,0.1,5H(O9A), 0.,5) 

C**RFHOVE  OUT-OF  RANGE-POINTS 

249  00  300  JJ»l,NUMOTS 
IF  (A(JJ).GE.SPLHTN)  GO  TO  25® 

A(JJ)*3.1254SPLHIN 

250  IF  (A<JJ) .LE.SPLMAX)  GO  TO  30n 
A{ JJI=SPLHAX-3.t25 
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besuvailadle  copy 


ion  ':onttijme’ 

FACH 

a (Mu^^Ti»i)  =«;JL  iiN 
a ( =S^L''-Trr 

XFO  (»1U*<  ■*^5  4-1)  =r 

XFKM  NU^PT 5» ?) =1. 
in  T = 1 , 

XFR(I)  =T  INC'FLOa-'  f T-i  ) 

CONTTNIJF 

f'ftLL  LIN'"(yF';,a,M'l»*P’';,l,1,  IT) 
VLOC  = YLO'-O.l‘= 

II=II+1 

IF  (IQ?a.*J<-.i)  r,n  791 

''ALL  D"' Al,  IM  (f  , M ■<,  S'M) 
in^Arin V -1 
in  ■'n  19''? 

•^qq  lONTiN'jr 

r,o  •'■0  ?oi 

399  CALL  OL0~  O.r', -3) 

r.o  10 

999  CO'^TINII'* 

CALL  PLO’E 
'5T0P 

CMH 


quqCQ’lTT'JF  GFTIN  (Y,  IPUNNC,  ICHMO,  A 
91-IENSIO'l  a (31) 

TF  (F.F9.  1 ) m ’9  E,p 
9FA9f9,  nOC)  O (L)  ,L=1 ,2) 

PEK^C.,  I 101  ) (A  (1.)  ,L  =■',l^) 

9FA  HI,  ID"!)  (A(L),L  = 1‘»,??) 
o*^af'('5,  ion-*)  (A  ID  ,L=?3,  ?“) 

10  00  FOP'^AT  (S’X,  ?FF.l) 

1001  Fno'^A''- ( lOFl.i  ) 

1003  format  I 1) 

GO  TO  -S 

50  PEAO  (5, ’0  00)  (A  (t.)  ,L  = 1,  1’) 

PEAO  (r,l0'’l)  (A  (L)  ,1  = 11 ,?':) 

RFAr  (•■,’001)  (A(L),L  = 21,?3) 

READ  (‘=,’00’)  GanAGF 

2000  FORMAT  (iex,13FG.l) 

2001  format  (1F6.1) 

2002  FORMAT(Ft.i) 

7F  WRITF  ((,,  ?0T) 

WRI’E(P,1  no>  (A  (t.),L=MMjNN» 
ion  F0RmaT(14  ,10E1’.P) 

200  FORMAT  (’XplOHOP 

OFTUPA) 

CRD 


fJN) 
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FUMCTIOM  nV‘‘!^AL 
A(3l) 

SUM=r . 0 

Tn  im  T=*^M,MN 

100  SUM  = S'JM4-(  to.''**  (AfT)/iC.n)) 

nvFc>CL  = i3  .n*ALnGn(su“) 

•’FTilPti 


SUR''CUTI"‘^  OGAtl^Mr  ,”“,NN) 

OIMrN<;Tn'l  A (3  1)  , WT  (??) 

OAT»i 

3— j—u«3f*3.'’f“t.9j'"0.3j'!.n>0.'3>l.njl.?jl.3>l.?jl.0»''«5y  — C.lf 

C-1. 1 

DO  mo  1=1, ?R 
IOC  A (T)  =A  ( T)  *Oi)AWT  ( T) 

WPI-'-f  (6,?f'n) 

WRITC  (e.,3n0)  ( A Cu)  ,l 

?n0  POOMAT  (■’X,lCHnoo  *••»**) 

30n  Fnr>MST  (IH  jlOFl'?.^) 
pFTDO‘1 

r,T) 

END  OF  PEO'I^’O 


STNFUT  SPLMTN=C0.,  FfNISH=?fi 
o 3LAPE,  lO'/C  •3)0  - 10  DFGOfcE  PITrH  <500 

7 

TOI  pni?  FOITS  044.0  0<-4.0  ')44.1  ni-4,')  044,0  Qa^.O  044.0  Q4-*.0  O48.O  043.5 
0*»9.7  051.0  O-jF.O  0f».7  '»''F.5  '’7R.0  064.?  07?. ■»  075.0  074.0 

0*5.0  0*''.7  076.?  076,"  P75.  7 07t,7  076.7  f'7(j.?  076.5  n76.5 

0*5.7  (J7  9.7  0 33.5  ^86. 7 "87.0  0e<0.7  054.*  044.0  044.  0 054.7 


*)t  WOM  C009*  049.  0 045.0  049.0  C‘.5.0  045.0  041.  0 043.  0 0 49.1  050.0  IfO.? 
155.5  0«3.7  0F6.7  0F4.7  073.0  031.2  073.5  004.2033.7  034.5 

■i*.«  "91.1  150.  1 090.7  C3.5  033.7  037.6  085.  0 034.7  0 34.7 

)•*.?  "9t.«-  154.2  094.'  0°9.6  095.0  0^2.5  049.0  049.0  1 03.7 

•I  CARP 

. -*»  • 
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BESr  AVyilMSiE  COPT 


v%\ 

0NE>TH1R0  OCTRVE  BRNO  CENTER  FREQUENCY  IN  HZ 
Figure  B4-1 : Plot  From  Somplo  NOISY2  Input  Data  Cords 
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I 


f 


APPENDIX  B-5 
PROGRAM  dBA 
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Ri8,STCSA,C'16  00  00  , 11*0,1050.  P720119  MCG.iEGO^ 
FTN, 

LGO. 


END  OF  RECORD 

PROGRAM  03A  ( I NPOT , OUTPUT,  PJNGH  , tape5  = IJPIJT  , TAPFt>=  OJ  Tr>U  T , T AP  c:7  r.ni  J 
ICH) 

INTEGER  INFO,FRQRNG 
REAL  JUN< 

DIMENSION  INF0(6),  03AAT(39),  08(2,40),  1MIKE(12) 

NAMELIST'INPUT/  IPUNCH 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

o 

c 

c 


DATA  0BA'(T/“80.0,“ij0.0,*60.*),*P0.0,“30.0,*70.**,*o3.^,“t>6.7  ,“?3.3,“ 
144.  7 ,*39.  4, *34.0, “30.2, “2b.  2,“22.5,*19.i,*lo.l,*i3.*,"10  . 3,*5.c-,*b 
2. &, *4. 8, *3. 2, *1.9, *0.6, 0.0, 0.6, 1.0,1. 2, 1.3, 1.2, 1.0, 0.3, *0.1, *1.1,* 
32. 5, -4. 3, *6. 6,-9. 3/ 

1 CONTINUE 

PROGRAM  08A  IS  USEO  TO  CALCULATE  30TH  03  ANO  D3A  NuiGHltJ  SOUNj 
LEVELS  FOR  ANY  SET  OF  OaTA  TAKEN  ON  THE  OUIET  PROFELUER  RUNS. 

THE  USER  MUST  SPECIPY  THt  NUM8ER  OF  <Pw  RUNS  PER  CoNFIGU  vhT 10  * A lO 
THE  NUMBER  OF  MICROPHONE^  INCLUDEJ  IN  THE  DATA.  T«0  JaTA  F)\-IhT: 

MAY  BE  read,  OEPENDIimG  ON  THE  uRI;jINAL  PUNCHFj  CArDs.  THu  PAIN 
OPTIONS  INCLUDE  CONSIDERATION  OF  ONLY  uESIREJ  MiCFOPH«)NES , A(U 
THEN  ONLY  OVER  THE  SPECIFIED  RANGE  OF  FREQUENCIES.  THt  Uj.Ek  )AY 
specify  THE  INTERNAL  PROGRAM  OPTION,  IPUNCH,  TO  Gl VE  JNlY  TnE 
PRINTOUT  ANU  NO  CARDS.  CAkOS  PROOUCEQ  3Y  THIS  PRuGRAn  NIcl  ^'>;-TAxN 
THE  09  ANO  DBA  LEVELS,  THE  i%UN  lOENTIFIJA  TIUN,  An  i HE  RANoE 
CONSIDERED  IN  THE  CALCULATIONS.  3ATA  CA>.J  Or^OERlNu  iS  GIVEN  TE;.o» 


SPECIFICATION  DATA  CAnJs 

NRJN  = NUMBER  OF  DIFFERENT  R°Ho/ I jNF IGURAT iON  iE 
NCHANL  = NUMBER  OF  MICkOPHUNE  CHANNELS/RPM  RUN  Is 
NMIKE  = TOTAL  NiJ.i3E\  UF  OESIREu  HiCROPHONEs  IE 
IMIKE  = MICROPHONE  IDCNTIFIcATIOh  NUMBER^  1015 
FRQRNG  = DATA  CAkO  FREQUENCY  RmNGE  CODE  Is 


1 . . 

3.15 

HZ  TO 

2C<  HZ 

c 

2 .. 

20  HZ  TO  20K  HZ 

c 

NN,  MM 

= KANGF  Tu 

BE 

CONSIDERED  T4 

CAL^'JaAT  I ONS 

215 

c 

NN  = 

S *■ 

INDEX 

OF  oThRTInG 

Fi^tOUCN'Y 

c 

MM  = 

3 ♦ 

TNOEX 

OF  final 

FREQUENCY 

c 

INDEX 

table 

c 

1 

HZ 

I 

HZ 

I 

HZ 

1 

HZ 

I 

HZ  I 

H7 

c 

1 

20 

6 

63 

11 

200 

lb 

63  0 

21 

200J  26 

c 

2 

23 

7 

MO 

12 

25  0 

17 

300 

22 

2530  27 

3)00 

c 

3 

32 

3 

IOC 

13 

315 

18 

1000 

23 

3150  2) 

1 0 0 0 0 

c 

4 

40 

9 

12^ 

14 

40  0 

19 

1250 

24 

■*JO0 

c 

5 

50 

10 

150 

15 

500 

20 

1630 

25 

3000 

C 

C 

c 

c 

<* 

KJ 

c 

c 

c 


PRESET  PARAMETERS  (NAMELIST) 

IPUNCH  s 0 ..NO  PUNCHtu  DATA  OtSIRtO 
(1)  ..  PUNCHED  DATA  OCSIREG 

• changes  TO  •^HtSE  Parameters  should  be  made  using  tme  jTanoai: 

NAMELlSr  FORMAT  (EXAMPLEI  ilNPUT  IPUNCH=3i> 


IPUNCHsl 
READ(5,  INPUT) 
RFAO  (3,19)  NRUN 
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READ  (5,19)  NCHANU 
READ  (5,19)  NMKE 

READ  (5,20)  (IMIKE(II) ,II=1,NHIKE) 

READ  (5,19)  FRQRNG 
READ  (5,21)  NN,MH 

MSTRT  = N-iIK£*-l 
□0  2 JJs.MSTRT,  12 

2 IHIKE(JJ)=99 
TLOBA=0.0 

3 READ  (5,22)  ( I NFO  ( 1 1 ) , 1 1 =1 , 9 ) 

IF  (IRFO( 1) .EQ.IOHRFSTART  ) GO  TO  1 
WRITE  (6,23)  (INFO( II) ,11=1,8) 

IF  (IPUNDH.EQ.O)  GO  TO  4 
WRITE  (7,22)  (INFO(II) ,11=1,8) 

4 IF  (IWFO(I) .EQ.IOH  ) STOP 

DO  18  L=1,NRUN 

ISTEP=1 

DO  Itt  I=1,NCHANL 
IF  (FRURNG.EQ. 1)  GO  TO  b 
DO  5 NMN=1,8 

5 D8(NMN|)=0.0 

READ  (5,25)  lORUN, lOMIKE , (09 (1 ) , N=9 , 16) 

GO  TO  7 

6 READ  (5,24)  I JRUr,  lO^'I  K£  , (0  3 (1 , N)  , N = 1 , 10  ) 

7 CONTINUE 

IF  (lOMiKE.EQ.  IMKE  (liTEP)  ) GO  TO  9 
DO  9 JJ=1 ,3 

8 READ  (5,26)  JUNK 
GO  TO  18 

9 ISTEP=ISTEPfl 

IF  (FRQRNG.EQ. 1)  30  TO  10 
READ  (5,27)  (Dd ( 1,N) ,N= 19, 26) 

READ  (5,27)  ( 06 ( 1 , N ) , N= 29, 36 ) 

READ  (5,27)  ( 3 B ( 1 , N ) , N=  3 9,  4 0 ) 

GO  TO  11 

10  READ  (5,27)  ( DB ( 1 ,N ) ,N= 1 1, 2 0 ) 

READ  (5,27)  ( OB ( 1 , N ) , N= 21 , 3 0 ) 

READ  (5,27)  ( 09 ( 1 , N) , N= 31 , 40 ) 

11  CONTINUE 

DO  12  J=NN,MM 

12  OB(2,J)=3D(l,J)i-D3AMT(J) 

DO  m K=l,2 

DO  13  J=NN,MM 

13  TL034=TL38At(10.**(05(K, J)/10.) ) 

D8(K,40)=10 (AtOGlO (TLU3A) ) 

14  TLOBA=0.0 

IF  (FKORNG.EQ.l)  GO  TO  16 

WRITE  (6,29)  09(1,40)  ,D6(2, 40)  ,13RUN,IJiIKE,NN,f1N 
GO  TO  16 

15  WRITE  (0,29)  08  ( 1, 4 0 ) , D3  (2 ,4  0)  , IORUN,  lOriKE  , NN,  N l 

16  IF  (IPUNCH.EQ. 0)  GO  TO  18 
IF  (FRQRNG.ED.l)  GO  '^O  17 

WRITE  (7,  28)  08(1,40) ,D9(2,hO) , ID, <UN,IOHIKE,NN, 41 
GO  TO  13 


17  WRITE  (7,  29)  08  ( 1 , 40)  ,D3  (2, 4 0)  , lURUN,  IPHIKE  , Ni),  NN 


IB  CONTINUE 
GO  TO  3 
C 

19  FORM/IT  (15) 

20  FORMAT  (1015) 

21  FORMAT  (215) 

22  FORMAT  (3A10) 

23  FORMAT  (1  HI , /// , 6 A 1 0 , / ) 

24  FURMAT  (oX,I2, 3X,I2,1X, lOF b.l) 

25  FORMAT  ( 3 X , 12 , 3X , 12 , 1 X , 1 JFfa . 1 ) 

26  FORMAT  (F6.1) 

27  FORMAT  (10F6.1) 

23  FORMAT  (2F10.2,30X,  BHIOI  7.  + x0,I2,3H  C0,I2,3H£  ,215) 

29  FORMAT  (2F10.2,30X,6HIOt  Ra,I2,3H  C0,12,3H?  ,215) 


END 

••  END  OF  RECORD 


JINPUT  $ 
6 
10 
4 

3 ‘ 

2 
9 


10 


36 


6-9LA0l£  bAStLINt 
IO«  •/.♦R006  COOli  01*7,7 

066.0  067.0  3a7.2  070. 

073.0  07P.0  031.0  0o6. 

044.0  094.0 

iOl  ■/.♦•R006  C002?  042. 5 
062.5  064.7  0S4.7  067. 

069.7  072.0  030.0  063. 

030.0  090.2 

IOl  •/.♦R006  C003i  042.5 

059.7  062.0  060.7  062. 

064.0  066.5  073.7  073. 

034. 0 035.0 

IOl  •/.♦R006  COOtJ  039.2 
C56.0  057.7  056.7  961. 
Obi. 2 064.2  071.7  076. 

029.0  062.0 


, 10  PITCH 

049,0  05-f.'’  058.2  377. J 090.5  065. 7 Obo.Q  079 
2 070.2  072.0  072.5  072.2  072.0  071.5 

0 062.7  363.2  065,7  079,7  04-t.O  0h4.J 

0‘*2.0  047.2  047.2  061.2  07h.5  0:^3. 0 0»3.5  Ool 
5 067,2  063.7  069.2  069,5  363.2  063.7 

7 079.7  065. C 084.7  077.2  039.7  03o.0 

041.3  0-+3.  2 043.2  0‘+9.0  0‘j7,7  OtS.j  052.7  061 
7 063.2  065.0  0o+.2  06h.2  063.2  063.2 

0 C7o.2  081.0  076.5  070.7  034.5  034.0 


♦ • ♦♦  ♦ ** 
.2  Of-  3.2 


.2  059.0 


.0  0 o . 0 


0h0.7  052.  0 043.  0 0-,3.2  057.7  J4j.2  J<»7. 
0 059.7  062.5  062.2  0o2.2  061.2  0o0.7 
7 073.0  077.0  074.5  066.0  029.5  029. 


05^.0  052,7 


200  SIMILAR  CARDS  REMOVED  TO  SHORTEN  LISTING 


IOl  Xi-KOll  C005S  O0O.2  r>56.0  0b4.  0 057.0  056.5  071.0  064.0  071.7 

078.5  091.2  082.0  093.2  091.0  091.0  090.2  090.5  090. j 009.5 

069.2  083.2  039.2  069.5  092.2  09h.7  097.7  100.3  067.5  057. 

054.0  104.5 

IOl  7.*^011  : 036»  064.2  059.2  055.0  061.  0 060.0  364.3  365.0  07  3.7 

378.5  094,0  032.2  091.7  099. 0 069.7  069.2  090.5  069.7  069.2 

089.2  087.7  036.2  388.5  090.5  093.7  096.7  u96.0  065.5  055.7 

053.0  103.5 

IOl  X«-R011  C007J  069.0  069.7  071.7  075.0  075.0  080.0  031.0  068.3 

093.7  098.2  092.2  105.0  099.7  104. » 101.5  103.5  102.5  101.0 

098.7  097,2  098.0  101.0  105.2  10o.7  103.0  110.7  073.7  U04.O 

064.0  115.0  * * 

IOl  •/.♦^Oll  0003$  066.5  065.0  0c5.  5 070.2  068.7  076.2  075.7  063.7 

089.7  097.2  090.0  103.7  102.2  101.7  101.7  101.7  103.2  098.7 

097.2  096.2  097.5  093.0  103,0  137.2  102.5  i09.0  079.0  064.0 

064.0  114.0 

lUI  7.»R011  G009S  056.0  058.0  058.0  064.2  0o2.2  0o9.7  070.9  079.7 

085.0  096.7  037.2  099.5  099.7  099.2  097.7  099.5  097.2  096.2 

095.2  095.7  096.2  095.0  099.2  IO4.O  099.2  103.7  374.0  056.0 

058.0  110.5 

IOl  7.»%011  COIOS  056.0  356.0  056.  0 060.2  058.7  069.3  073.2  074.2 

080.7  095.2  034.  2 395.2  0<>2.  7 094.2  094.  2 093.5  092.2  092.0 

091.7  091.7  091.5  091.0  093.0  093.2  094.0  098.0  0b7.0  056.0 

096.0  105.2 

— blank  card  


0o9.7  0,  7.2 


092 .V  0,  3,5 


101.  J 0''-3.5 


100.2  090.7 


098.5  0o3.0 


091.7  OjI.O 
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♦♦  6-^LftjE 

9AStLlNt,  10  Otb^tE  PITCH  ***• 

. * < 

* 

85.tJ0 

85.66 

lot 

•/.•*•  XO 

o 

CO  3.* 

ib 

82. &3 

82.73 

131 

Z*kO 

6 

CO  4i 

93.13 

93.22 

lOI 

7.  + R0 

3 

C 3 9^ 

y 

i :i 

83.23 

88.30 

10  t 

•/.♦HO 

O 

COlO? 

y 

5b 

91.09 

90.39 

lot 

•/.♦«0 

7 

CO  3* 

y 

7 - 

87.25 

86. 68 

lot 

■'i  + KO 

7 

CO 

.j 

3t_» 

99.70 

99.0  0 

nt 

■/.  ♦XO 

7 

C 0 

9<».95 

93.99 

lot 

■/.♦kO 

7 

COIOS 

j ? 

9^.87 

94.26 

lot 

•/.♦RO 

8 

CJ  35 

90.71 

90.16 

lot 

■/.♦kO 

i 

w3  ti 

y 

7 . 

103.53 

102.88 

lot 

•/.♦ro 

i 

6 J 91 

2fc 

97.95 

97.39 

lot 

•/.♦xo 

3 

cnos 

y 

i-,9 

9ii.67 

95.91 

10  t 

•/.♦RO 

9 

CO  3i 

y 

3c 

93.02 

92.37 

lot 

•/.♦p,o 

9 

CO  ■»£ 

3o 

107.38 

106.47 

lot 

•/.♦KO 

CO  3J 

t ^ 

101.64 

100.61 

lOt 

Z+P.O 

9 

k/  0 1 0 ^ 

y 

: c. 

100.41 

99.24 

nt 

•/.♦XOlO 

CO  3i 

y 

96.87 

95.84 

lot 

•/♦ROlO 

CO  -»$ 

i •. 

109.93 

109. C4 

lot 

•/.♦xOlO 

00  95 

y 

r> 

104.42 

103.61 

lot 

/.♦R310 

■^0105 

y 

103.65 

101.55 

lot 

•/.♦xJll 

CO  3> 

y 

j z. 

99.5  0 

97.82 

lot 

CJ  -.J 

• / 

APPENDIX  B-6 
PROGRAM  NEWCITY 
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Rt13,STCSA,C'155ij00,T3C,IO60.  3720119  MCGREGOR  »5h21 
FTN,R=3, 

HAP, PART . 

ATTACH,CCAUX,C-AOX, I0=X654321. 

ATTACH,  CC6600 , 1 D = Xb5<+ 521. 

LI3RARV,CGAUX,:C6600. 

LGO. 

“ END  OF  RECORD 

PROGRAH  NEHCiTY  (INPUT, OUTPUT , TAPE5=INPJT, TAPE6=0JTPuT , PLOT) 

C 

C PROGRAH  NENCITY  IS  DESIGNED  TO  PLOT  SPL  »/S.  HORStPOHER  AND  TnRUST 

C FOR  ANY  PROPELLER  RUNS.  THE  USER  HAY  SPECIFY  THE  UiE  OF  ONLY 

C PART  OF  THE  DATA  FOR  PLOTTING  AND  HAY  RECEIVE  THE  DATA  IN  EITHER  D3 

C OR  D3A.  CURVES  MAY  Bt  EITHER  SPLINE  FITTED  OR  UP  iO  jIXTH  OkOlR 

C POLYNOHIAL  FITTED.  NAHELIbTED  PARAMETERS  ARE  LISTED  3ElOR  ALUNG 

C WITH  instructions  FOP  THEIR  CHANGE. 

C 

C *•  DATA  TO  3E  READ  IN 

C SPECIFICATION  OF  NUH9ER  OF  RPM  VALUES  IN  RUN 
C NRPM  = RPM  VALUES  PER  RUN 

C 

C SPECIFICAT ION  OF  HICROPHONES  IN  DATA 
C NMIXE=NUMBER  OF  MIKES 

C I0HIKE=M1CR0PH0NE  IDENTIFICATION  NUMBERS 

C PRESET  PARAMETERS  (NAMELIST) 

C 

C N(iJRV£  = (0)  ..  SPLlNt-FIT 

C 1 THRU  6 ..  NTH  ORDER  POLYNOMIAL  FIT 

A 

W 

C IOPlOT  = 1 ..  HORSEPOWER  ONLY 

C 2 ..  THRUST  ONLY 

C (3)  ..  BOTH 

C 

C LEGEND  : (1)  ..  LEGEND  BLOCK  PRINTED 

C 0 ..  BLOCK  OMITTED 


C 

C 

c 

c 

c 

c 

c 

A 

w 

c 

c 

c 

c 

c 

c 

c 

c 


lORPM 


(1)  ..  ALL  available  DATA  USED  IN  PLOTS 
2 ..  3 .RPM  VALUES  PLOTTED  REGARDLESS  OF  DATA  INPUT 


lUOBA  s (0) 
1 


..  OUTPUT  IN  DB 
..  OUTPUT  IN  DBA 


YHIN  : (40.)  ..  HININJH  06  OR  DBA  FUR  PLOT 

YMAX  : (120.)  ..  MAXIMUM  06  OR  DBA  FOR  PLOT 

XMINT  « (0.)  ..  MINIMUM  THRUST  FOR  PLOT 

XMAXT  * (1625.)  ..  MAXIMUM  THRUST  FOR  PLOT 

XHINH  > (0.)  ..  MINIMUM  HORSEPOWER  FOR  PLOT 

XMAXH  « (260.)  ..  MAXIMUM  HORSEPOWER  FOR  PLOT 


» CHANGES  TO  THESE 
namelist  FORMAT 


PARhMETERS  should  be  made  using  the  STANDARD 
(EXaMPLEI  SINPUT  NCURVE*2,  XMINT*40.* 


INTEGER  MEAD 
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DIMENSION  HEA0(5),  X(2,3,5),  Y(4,3|5),  XMIKEdO),  LMlKECt),  llTtEX 
1(2, 2>,  ritUEV(^),  rPL8Ll(2),  BETACi),  XSIEP(2>,  XDAJ(2),  yjAT(7), 
aiOMIKEdO),  :x(in3),  CY(i03),  W(5),  WORK(IOO),  CO£F(7),  ORJ(b), 
3XM1N(2) 

NAMELIST/  INPJT/  NCUR VE , lOPL  JT  , LEGEND , IORP.1 , 10 J d A,  Y .IlN,  Y r1  AX , XMI  NT  , 
IXMINH, XMAXr ,XMAXH 

DATA  TIT-EY/IOHSOJNO  LE(/E,10HLS  IN  OB  R,iOHE  0.0G02  H,10HICkO3Ar 
1 /,  (7ITlEX(1, I) ,1  = 1,2) /liH  HORSE, lOH  THRUST/, (TITlEX ( 2 , 1)  , i 

2=1,2)/10HPOWEP  ,10H  (.dS)  / , TPlBl  1/ lOHMI  KE  NUr13E,10HK 

3 /,  0R0/5HlST,3H2N0,3H3R0,3H^TH,3hSiH,3HbTH/ 

iO\/E  PEN  TO  NEW  ORIGIN 
CALL  Pl.0T  (0.,1.,-3) 

NCJRVE  = 0 $ IDPL0T  = 3 i LEGZNU=1  $ IORPM=l  $ I0D3A=d  f YilN=s*0. 

YMAX  = 120.  £ Xf1INT=0.  I XMI.mH  = 0.  £ XMAXT  = lb25.  ! XMAXH=26tj. 

1 READ  (5, INPUT) 

READ  (5,23)  NRPN 
READ  (5,28)  NHIKE 

READ  (5,29)  ( lOHlKF ( JJ) , JJ=1,NMIKE) 

NLINE=1 

IF  (NDJRVE.NE.  0)  NLINE  = -1 
DO  2 JJ=1,NMIKE 

2 XMIKE( J J) =FLOAT (iDMIKE ( JJ) ) 
rSTEP= (YRAX-YMIN)/3.0 
XSTEP(l) = (XNAXH-XMINH) /o  .5 
XSTEP(2)  = (XMA  XT-XMINT) /o.5 
XNIN(l) sXHINR 

XHIN(2) =XrtINT 

••  read  in  HuRSEPOWER  and  THRUSi  VALUES 

3 CALL  Reader  (head, beta, x,y,nkpm,nnike,ijrpm,NjIZ,iod3a) 

SJMMA=0.a 

LINO=0 

if  (HEA0( 1) .EQ.lOH  ) GO  TO  27 

IF  (HEAD( 1) .EQ. ICHRESTaRT  ) GO  TO  1 

IF  (Ioplot.ne.i.ano.iuplot.ne.2)  go  to  a 

LIN0=1 
L=I0Pl0T 
go  to  5 

*♦  continue 

u = l 

5 CONTINUE 
b DO  26  I=1,NMIKE 
DRAW  AXIS  AND  LABELS 
IF  (1009A.EQ. 1)  GO  TO  7 

call  axis  (0.,0.,TITLEY,fa,3.,90.,TMIN,Y3TEP) 

GO  TO  6 

7 CALL  AXIS  ( 0.  , 0. , 38HSDUN0  lEVELS  IN  DBA  RE  0.00J2  .ilCROBAR,  3 3 , o . , 9 

10«,hQ«,10*) 

8 CALk  AXIS  (0. , 0.,TITlEX(1,l) ,-20,6.i,0.,XMIN (L) ,XSTEP(L) ) 
call  SYMBOL  ( 2. 25 , 8 . L , 0 . 15 , TFl BLl , 0 . , 20 ) 

CALL  NUMBER  ( 4 . 05 , 8 . A , 0. 15 , XMlKE ( 1) , 0 . , -1) 

CALL  SYHBUw  ( 0 . 85 , 3 . 25, 0 . 1 , HEAD , 0 . , 50 ) 

IP*0 

DRAW  CONSTANT  BLADE  mNGLE  CURVES 

9 IJ=1 

jPOiNTsNLlNE 
IN0  = 0 


I 


I 


c 


c 


J=1 

10  30MTIMUE 

00  11  K=1,NSIZ 
YOAT (<) =Y 

11  XOAT  (<)  =X  (L,  J,  K) 

SKIP  5 OEGKEE  DUHMY  UATA  LINES 
TOT=TDTAl (YJAT,NSIZ) 

:UTOFF  = *tO.O*FLOAT(NSIZ) 

IF  (TOT.lE.CJTOFFI  go  to  23 
JK=J 


EBi  Av/am  m 


(=  tO.  Dtj/U3A) 


30  TO  15 

••  DKAW  CONSTANT  RPN  CURVES 
12  IJ=2 

JPCINT=0 


IND  = 1 
< = 1 


13  CONTINUE 

□0  14  JJ=1,3 
YOAT ( J J) =Y (I, JJ,K) 
lit  K0AT(JJ)  = X(L,  JJ,K) 

C SKIP  any  OJNiY  LINES 
T0T=T0T4. { YDAT,3) 

IF  (TOT. LE. 120.0)  GO  TO  23 

JUNP=3 

JK=K 


50  TO  lo 

15  ISYN=J-1 
JU1P=NSIZ 

16  XOAT(JUl1P  + l)  =XMIN(,_) 

XDAT (JJNPf 2) =XSTE?(L) 

YOAT ( JJN^^l) = YMIN 
YOAT (JJN3f2)=YSTEP 

C ••  REMOVE  5 jESKEE  OUMMY  DATA  POINTS 
IF  (YOAT (2)  .jT. 40.0)  GO  TO  Id 
MOVE=JJMPfl 
00  17  KI=2,M.<VE 
XOAT  (KI)  = X0AT  (KIH) 

17  Y0AT(KI)=Y0AT(Kin) 

XDAT(MUV£*1)=0.0 
YOAT  (MOVE^DsO.O 
JJMP=JJMo-l 

C ORAM  FITTED  CURVE  AND/OR  jATA  POINTS 

18  IF( JPOINT.EQ.O.ANO.IJ.LO.2)  GO  TO  19 

SALE  FLINE  (XOAT, YOAT , -JUMP, ItJPOINTf ISYH) 

YLOC=0 .6- (.15*FL0AT (I5YMI) 

IF  (LEGENJ.Ea.O.OR.IJ.EQ.Z)  GO  TO  19 

C DRAM  legend  3LOCK  FOR  PARTICULAR  LINE 

CALL  5YM30L  ( 4 . 20 , YLOC , 0 . 1, 22H--  DEG.  3LAUE  ANGLE, 0. ,22) 

CALL  5YM30L  ( 4 . 0 , ( YL OC t . 05)  , U . 1 , ISYM , 0 . , - 1) 

CALL  NJM3ER  ( 4. 5 , YLOC , 0. 1, 3ET A ( J<) , J . , -1 ) 

19  IF  (NCJRVE.EU. 0)  GO  TO  23 

C CALCULATE  AND  DRkW  LEAST  SQUARES  FITTtO  CURVE 

00  20  11=1, JUMP 
CXI  ID  =X0AT  (II) 

20  CY( II) »Y0AT (II) 

STEP=(CX(JJM?) -CX(1) ) /lOP. 
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MCP  = NCJ^\/E 
ICPk=N0U<(/E  + 1 
IF(  JUMP._T.I3PR) 

PRINT  *,'*pi\iS  "fJ,"  through  PLSC'=’“ 

CALL  PLS:F  (LX,CY,-l,jJMF,N:F,NMAX,COLF,0,XD,<O,»^ORK,iLt<) 

PRINT  *,  “nCP  = ",NCP,”  EKRUR  NUM3ER  = ’MER,"  JMAX  = 'SN.ImX 

PRINT  *,‘*:OErS  = ",  (CDEF  (<K)  ,<K=l,  7) 

i_UPEND=>HAX+l 

DO  22  JJ=1,1G1 

3X ( J J*1 ) = CX ( JJ) +STEP 

DO  21  <<=1,LJPEND 

21  5UMMA  = S JMMA  + (C0EF(KK)*(CX(JJ)  ♦♦(«-!)  ) ) 

CY( J J) =S JMMA 

22  SUMMA=0.n 

CX  ( 102) =XMIN (L) 

CX( 103) =XSTEP(L) 

CY (102) =YrtIN 

CY(103) =YS1 EP 

PRINT  • , ( CX (KK) ,KK=1,1C3) 

PRINT  *,  (CY(KK)  ,i<K=l,103.) 

CALL  LINE  (CX,CY, 101,l,u,iSYM) 

PRINT  »,"LINE  DRAWING  COMPLETED" 

3RINT  *,  " " 

IF(NCURl/E.EQ.O.OR.IP.tC.l)  GO  TO  23 

CALL  SYMBOL  (2  . o , 0 . 1 , 0 . 1 0 , 37  H*  ♦ ORDER  POLYNu-IIAl  FxTTEj  CUR7E3 

1,0. , 37) 

CAlL  SYM30L  ( 2 . 9 , 0 . 1 , 0 . 1 , OR D ( NCP ) , 0 . , 3 ) 

IP=1 

23  CONTINUE 

IF  (IND.EO.l)  Gj  to  24 
J=J*1 

IF  (J.LE.3)  GO  TO  10 
30  TO  25 

24  CONTINUE 
< = <♦1 

IF  (K.LE.NSIZ)  GO  TO  13 

25  CONTINUE 

IF  (lU.EJ.l)  GO  TO  12 
C MOVE  »EN  TO  NEW  ORIGIN 

2b  CALL  PuOr  (0.5,0., -3) 

IF  (LlwO.EO.l)  GO  TO  3 
L=L*-1 
LINO=l 
SO  TO  5 
27  continue 
CALL  PLOTE 
STOP 


28  FORMAT  (15) 

29  FORMAT  (1015) 
END 
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SUB^OJTIME  t^EADEK  (H£AO,t3ETA,X,Y,i\(ROM,NMIi^E,  IOKP.“1»  NS±Z  , I003A) 
DIHENSIOM  HEA0(5),  X(2,3,j),  Y(h,3,5),  9ETA(3) 

INTEGER  3EA0 
READ  (5,14)  HEAD 

IF  (HEA3( 1) .EQ.lOH  . OR .HEA D ( 1 ) . EQ. 1 OHREiT ART  ) GO  TO  13 

IF  (IJRPM.E3.1)  GO  TO  9 
IF  (IORPM.EQ.2.AND.NKPM.EQ.3)  GO  TO  9 
C READ  3 POINTS  IF  3 ARE  AVaIlABLE 
DO  2 1=1,3 
READ  (5,15)  SETAd) 

DO  1 J=l,2 

READ  (5,1b)  RPM,X(1,I,J) ,X(2,I, J) 

1 READ  (5,17)  TRASH 

2 READ  (3,16)  RPM , X ( 1 , I , 3 ) , X ( 2 , I , 3) 

IF  (I0D9A .EQ.l)  GO  TO  5 

C *•  READ  03  DATA 

DO  4 I = 1,NMIKE 
DO  4 J=l,3 
DO  3 <=1,2 
READ  (5,17)  Y(I,J,<) 

3 READ  (5,17)  TRASH 

4 READ  (5,17)  Y(I,J,3) 

GO  TO  8 

C ••  READ  D3A  DATA 

5 DO  7 I = 1,NMKE 
DO  7 J=l,3 

DO  6 <=1,2 

read  (5,15)  Y(I,J,<) 

5 READ  (3,18)  TRASH 

7 READ  (5,18)  Y(I,J,3) 

8 N5IZ=3 
SO  TO  13 

C *•  READ  A.L  AtfAI.ABUE  DATA 

9 DO  10  1=1,3 

READ  (5,15)  3ETA(I) 

00  10  J=1,NRPM 

10  READ  (5,16)  RPM,X(1, I, J) ,X (2,1, J) 

DO  12  I = 1,NMKE 

DO  12  J=l ,3 
00  12  <=i,NRPM 
IF  (10D3A.EQ.1)  GO  To  11 
C REAj  D3  DATA 

READ  (5,17)  Y(I,J,<» 

SO  TO  12 

C ••  READ  03A  DATA 

11  READ  (5,18)  Y(I,J,<) 

12  DONTINJE 
NSIZ=NRPN 

BEl  AVAIWDIE  COPY 


13  RETURN 

14  FORMAT  (8A10) 

15  FORMAT  (3uX,F10.3) 

16  FORMAT  (3F10.3) 

17  FORMAT  (rlO.3) 

18  FORMAT  (10X,F10.3) 
END 
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